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Abstract 

The Stokes multipliers in the matrix models are invariants in the string-theory 
moduli space and related to the D-instanton chemical potentials. They not only 
represent non-perturbative information but also play an important role in con- 
necting various perturbative string theories in the moduli space. They are a key 
concept to the non-perturbative completion of string theory and also expected to 
imply some remnant of strong coupling dynamics in M theory. In this paper, we in- 
vestigate the non-perturbative completion problem consisting of two constraints on 
the Stokes multipliers. As the first constraint, Stokes phenomena which realize the 
multi-cut geometry are studied in the symmetric critical points of the multi-cut 
two-matrix models. Sequence of solutions to the constraints are obtained in general 
fc-cut critical points. A discrete set of solutions and a continuum set of solutions 
are explicitly shown, and they can be classified by several constrained configura- 
tions of the Young diagram. As the second constraint, we discuss non-perturbative 
stability of backgrounds in terms of the Riemann-Hilbert problem. In particular, 
our procedure in the 2-cut (1,2) case (pure-supergravity case) completely fixes the 
D-instanton chemical potentials and results in the Hastings-McLeod solution to the 
Painleve II equation. It is also stressed that the Riemann-Hilbert approach realizes 
an off-shell background independent formulation of non-critical string theory. 
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1 Introduction 



Non-critical string theory [T] has provided interesting theoretical laboratories which un- 
cover various intriguing features about string theory. This string theory is known as solv- 
able system not only in the perturbative world-sheet formulation, Liouville theory p]-[TO] , 
but also in the non-perturbative matrix-model formulation pn432] . Recently, among 
various kinds of matrix models, the multi-cut matrix models |33] have turned out to be 
a fruitful system. The first discovery was on the two-cut matrix models [3H439] . which 
were found to describe type superstring theory [10^H2] . Furthermore, the multi-cut two- 
matrix models were generally found to have a correspondence with the so-called fractional 
superstring theory [l3] and also with non-critical M theory as its strong-coupling dual 
theory [H], which realizes the philosophy proposed in the Hofava-Keeler non-critical M 
theory [i5] . 

Quantitative analyses of critical points and perturbative amplitudes in the multi- 
cut two-matrix models have been carried out in [mHB]. The main observables used 
there are macroscopic loop amplitudes (or resolvent) p^] - [T5lfT71fT8l[29l[30lH71 - [5T] which 
provide the information of spectral curves, the classical spacetime of this string theory 
[32l[52l[53]. A concrete expression for spectral curve is important because it provides 
relevant information for reproducing all order perturbative amplitudes in the multi-cut 
two- matrix models by the method of topological recursions [M]. 

The main theme in this paper is, on the other hand, about non-perturbative aspects 
of the multi-cut two-matrix models. Non-perturbative aspects in matrix models have also 
been studied extensively pij [271 [281 1311 [52|l51ti53t[55M80jFI The main concern is about 
non-perturbative contributions to the matrix-model free energy J^{C] fistr) on the large 
spectral curve 

oo 

J^{C; fi'str) ^ V fisto~^ 

•^niC) ~l~ J' non-perturb. (C; fi'str), fi'str -^0. (1.1) 

asym ' ' 
n=0 

Here J>i(C) is the genus-n perturbative free energy on the spectral curve C, and the 
information of the matrix-model potentials (so-called KP flows {tn}ngz) is implicitly 
included in the spectral curve: 

C = C{{t^}ne^), {t„}„ez e A^ii~^ c C-. (1.2) 

Here A^ltHng """^"^ stands for the non-normalizable string-theory moduli space mE The 
first quantitative implication was given in the early 90 's and is about the strength of string 
non-perturbative corrections which are of order (9(e~^/^^*'') quantities [82], i.e. open-string 
(D-brane) degree of freedom [83] : 

-^non-perturb. (C; fi'str) = ^/ exp S^^ti^'i 9str) ■ (1-3) 

r '- fi'str -I 



^See f80^ for a nice review of these recent progress. 

^We carefully put "asym" below the equation in order to emphasize that they are equal only in the 
asymptotic sense. 

•^The normalizable string-theory moduli space Mitring ^ known as the space of filling fraction f76] 
which parametrizes the on-shell string backgrounds. The off-shell backgrounds are defined in Section [5] 
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Here J is a set of indices which labels multi-instanton sectors, / = {ii, i2, - ' ' }y 



9str) = J2 ^inst(C) + ^(9str). (1-4) 

i&I={ii,i2,--- } 

Each primitive instanton action sl^^^{C) {i = 1,2, ■ ■ ■ , iVinst), is shown to correspond to a 
singular point of the spectral curve C [i2 l[56l - [60ll62ll69] and is identified with the ZZ-brane 
disk amplitudes in Liouville theory [3191110] • It is worth mentioning that these instanton 
corrections including higher order Qstr corrections S^^ti^] 9stv) are generally expressed 
as theta functions on the spectral curve [55[[73] and important in order to make the 
free energy J^{C; Qstr) modular invariant under modular transformations of the spectral 
curve C and also to be background independent in the normalizable string-theory moduli 
space A^strin™ '' the filling fractions) [73l[76]. The constant 6i is called D- instanton 
chemical potential (or fugacity). These constants are understood as integration constants 
of corresponding string equations [28], that is, 

fi = o, mez, Mnez e M^:zr"'-\ (1-5) 

for the flows in the non- normalizable moduli space A^string"^™™ shown [32] 
that the only A^inst (i-e. the number of primitive instantons) chemical potentials 6i {i = 
1, 2, ■ ■ ■ , A'inst) are independent among all the chemical potentials 6*/. 

Although various aspects of matrix models have been understood well so far, there 
still remains an important issue regarding the D-instanton chemical potentials. This 
is also known as non-perturbative ambiguities of string theory. Therefore, what is the 
physical requirement to determine the D-instanton chemical potentials? Although the 
actual matrix models should employ some particular universal values [01], they seem 
to be totally free parameters at least within continuum formulations based on string 
(or loop) equations. This point has been studied in the bosonic minimal/2D string 
theories [6TlE5l|66l[68] , in the type (1, 2) superstring theory [63], in the collective string 
field theory [63], in the free-fermion formulation [511 [67] and in the topological string 
interpretations [76]. In this paper, we address this issue by solving non-perturbative 
completion problem within a continuous formulation for the critical points of the multi-cut 
two-matrix models. In practice, we pick up physically acceptable D-instanton chemical 
potentials which realize physically reasonable behaviors in the non-perturbative regime 
(7str oo. Our solutions are based on two physical requirements: One is multi-cut 
boundary condition (in Section Hj) and the other is non-perturbative stability of perturbative 
backgrounds (in Section E]). 

The first requirement, the multi-cut boundary condition, is a non-perturbative con- 
straint on the Baker- Akhiezer function system in these multi-cut critical points: 

gstrg^^it- = Q(t; C) ^{t- 0, ^?str^^(t; C) = V{t; () ^{t- Q, (1.6) 

where the equation system here is expressed as an ordinary differential equation in ( and 
its isomonodromy deformation system in tE Note that the Lax operators in Eq. (11.61) 



The parameter t is one of the parameters in the non-normahzable moduli space A^st^Tng which 
is usuahy a couphng of the most relevant operator or the world-sheet cosmological constant. 
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in the fc-cut critical points are k x k matrix- valued operators (BS]. The idea of the first 
constraint is motivated by the non-perturbative relationship between the Baker- Akhiezer 
functions and cuts in the resolvent curves. This kind of relationship is discussed in terms 
of Airy function |52] . Specifically, the asymptotic expansion of the Airy function around 
the cut {( — i- — oo) is expressed a^ 



Ai(t;C) ^ f ,/""' J ' e-^('^+*)'''+^e5^(^+*)'''" +■■•, (1.7) 



a^m V(C + t)l/2 

where the relation to the resolvent (or macroscopic loop) operator 7^(C) 120] is roughly 
expressed as 

with the expectation value (■ ■ ■) which is taken with respect to the Hermitian one-matrix 
model of a matrix M. From this expression, one observes that the cut in the negative 
axes {( < —t) appears as a line where a competition between the exponents e ^sstr^^ 
(i.e. along the Stokes lines) happens. Therefore, we interpret this as a non-perturbative 
definition of the resolvent cuts. This consideration turns out to be important in the 
fractional- superstring critical points of the multi-cut two-matrix models since most 
of the cuts in these critical points are created by this procedure and cannot be read 
from the algebraic equations of the resolvent spectral curve. However, as we will see in 
Section HI this procedure do not necessarily create the necessary and sufficient k cuts on 
the resolvent curve, even though the fc-cut Baker- Akhiezer function Eq. ( 11. 6p is obtained 
from the assumption that the critical points have k cuts around ( ^ oo. In view of this, 
we need to impose a physical constraint so that the resolvent curves in the k-cut critical 
points should have k cuts around ^ — )■ oo. This constraint is expressed in terms of Stokes 
multipliers for the possible Stokes phenomena in this system. 

The second requirement, the non-perturbative stability of perturbative backgrounds, 
is imposed in the other formulation which is closely related to the Baker- Akhierzer func- 
tion system: the so-called the Riemann-Hilbert (or inverse monodromy) approach [Ml - 186] 
|23j . A brief flowchart of this approach is shown in Fig. [H Details are given in Section |5l 
but in order to show how the Riemann-Hilbert approach works in resolving the issue, we 
here show the leading expression of the free energy (more precisely the two-point function 
of cosmological constant t) in the two-cut (1,2) case: 

= [mY^ m = / |^e^"'(*^^)-'^'(*^^) + ■ ■ ■ . (1.9) 

The parameter Sn,2,i is a Stokes multiplier of the Baker- Akhierzer function system of the 
corresponding integrable system and the contour /C„ is an anti-Stokes line corresponding 
to the Stokes multiplier s„^2,i- As one can suspect from the expression, the Riemann- 
Hilbert approach is directly related to the study of Stokes phenomena at C — ^ oo in the 
ordinary differential equation of the Baker- Akhierzer system. 

In this expression, the function g^^\t; Q is an arbitrary function hut should be properly 
chosen so that the integrals other than the "leading" expression shown in Eq. (11.91) are 



^The asymptotic expansion of Airy function is reviewed in Appendix \X\ 
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Figure 1: The Riemann Hilbert approach and the D-instanton chemical potentials (or fugacities) 



negligible [86]. From the matrix- model viewpoints (to be discussed in Section [5]), this 
function can be interpreted as an off-shell string background geometry of string theory. 
Therefore, if one chooses g^^\t] Q as a macroscopic loop amplitude realized in the large N 
limit of the matrix models, then the leading integral (11 .Qp becomes a similar expression 
to the mean field expression for a single eigenvalue of the matrix integral which appears 
in various studies in literature [28l|31],[32],|60l|6Tl[63j|^ Therefore, the Stokes multipliers 
Sn,2,i in Eq- (11 -Sp are directly identified as the D-instanton chemical potentials in the 
semi-classical saddle-point analysis. That is, the first constraint is directly related to 
the constraint on the D-instanton chemical potentials. Furthermore, since the Riemann- 
Hilbert integral, Eq. (11.90 . provides the complete integration representation based on 
the reference string background g^-'\t;C), we can discuss non-perturbative stability of 
the background g^-'\t;C), especially for the background which is obtained as large N 
limit of the matrix models. This consideration for the stability is also expressed as a 
constraint on the Stokes multipliers and therefore the D-instanton chemical potentials. 
Originally, the mean field analyses include ambiguity of choice of contour and weight 
of these contours [28] and this fact becomes a cause of the ambiguity about the D- 
instanton chemical potentials in continuum loop-equation systems. In the Riemann- 
Hilbert approach, however, these degrees of freedom are identified as anti-Stokes lines /C„ 
and Stokes multipliers Sn,2,i, and they are tightly related to each other. As a consequence, 
the physical section of the D-instanton chemical potentials is obtained in the name of non- 
perturbative completion. This viewpoint is important in non-critical string theory because 
non-critical strings are sometimes defined as the large N (i.e. perturbative) expansion 
of unstable matrix-model critical points (e.g. (2, 3) bosonic minimal string theory) and 
therefore the matrix-model description does not necessary guarantee non-perturbative 
completion of string theoryQ 

As we will see in the coming sections, the above procedures completely determine the 
D-instanton chemical potentials in the two-cut (1, 2) critical points and results in the 
Hastings-McLeod solution [89] to the Painleve II equation (in Section ISTTl) . Actually it is 
known that this is the unique solution which realizes the two phases of the two-cut (1,2) 

^It is interesting that the Riemann-Hilbert expression gives a similar expression to the D-instanton 
operators obtained in the free-fermion formulation |31ll32j . 

^Early investigations of non-perturbative complete string theories are found in [TOIIMI 1571155] . 
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critical point of the two-cut matrix modeljj and therefore the Hastings-McLeod solution 
is suitable for this critical point. An advantage of our work is the discovery of the actual 
physical requirements to obtain the correct solutions to the non-perturbative completion 
which are also applicable to the critical points with an arbitrary number of cuts. 

The next main developments shown in this paper is, therefore, an extension of our 
procedure to the general multi-cut cases (which even reaches to oo-cut!). In particular, 
general structures of Stokes multipliers in the k x k isomonodromy systems are investi- 
gated in Section |3l and a new way to identify non-trivial Stokes multipliers is proposed 
(Theorem!!]) with terminology of the profile of dominant exponents. Furthermore, explicit 
solutions are obtained with help of the physical constraints, i.e. the multi-cut boundary 
conditions (Theorem [7] and [H] in Section [4.3p . In this sense, our solutions provide the 
multi-cut generalization of the Hastings-McLeod solution. Interestingly, we found that 
these solutions are labeled by constrained Young diagrams (Proposition [3] in Section H73|) . 
This result implies that there is a quite rich world beyond the non-perturbative horizon, 
and that the multi-cut matrix models provide fruitful fields for a quantitative study of 
these issues. 

Organization of this paper is as follows: In Section [21 after summarizing the asymp- 
totic expansion of the ODE system in the multi-cut critical points, the general facts 
about Stokes phenomenon in ordinary differential equations are reviewed. As a warming 
up, the case of the two-cut (1, 2) critical point is also shown. In Section [3l Stokes phe- 
nomena in the multi-cut critical points are studied. In particular, a systematic way of 
reading the Stokes multipliers in general cases is developed. In Section HI the multi-cut 
boundary condition is proposed. In Section 14.31 the discrete and continuum solutions are 
shown. In Section [5l the non-perturbative stability condition is studied in terms of the 
Riemann-Hilbert problem. Section E] is devoted to conclusion and discussion. 

Context of Appendices is: Appendix |X] is about the Stokes phenomenon of Airy 
function (a review of [52]). Appendix [B] is about calculation of Lax operators. Appendix 
[Clis about supplements to Theorem [5] and Theorem [6] with some examples of the multi-cut 
boundary-condition recursive equations. Appendix [D] is about derivation of continuum 
solutions. Appendix [E] is about calculation of the 3-cut (1,1) critical points and Appendix 
[Fjis about calculation of 4-cut (1, 1) critical points. 

2 Stokes phenomena in the ODE systems 

Before we devote ourselves into the multi-cut systems, here we first review some gen- 
eral facts about Stokes phenomenon in ordinary differential equation systems, then we 
summarize the well-studied two-cut (1,2) case. This two-cut system has been exten- 

^It was shown by Hastings-McLeod [5^ that their solution is a unique solution to the Painleve II 
equation, Eq. (|2.43|) . which realizes the following asymptotic behaviors of f{t) on the two sides of infinity 
t —i' ±oo: 

lm^f{t) + 2tfit) = 0: /(i^oo)^O, f{t^-^)^V~t, (1.10) 

which is the same behavior as the two-cut (1,2) critical point of the two-cut matrix model discussed 
in [42]. For some mathematical derivation of this solution in the two-cut matrix models, see also [90] 
which has been studied within the Riemann-Hilbert problem. 
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sively studied not only in physical context [5^HS5| I ^ [55 | I5U] but also in mathematical 
context [85 | [86 | [89 | [921495] ■ since it is related to the Hastings-McLeod solution |89] of the 
Painleve II system. For more comprehensive and rigorous reviews and references on 
the isomonodromy deformations, Stokes phenomenon and inverse monodromy problems, 
see [91]. We also note that the idea of ismonodromy deformation was introduced in 
non-critical string theory by [23] . 



2.1 The ODE system and asymptotic expansions 



It was first proposed in [69| that the multi-cut matrix models are controlled by multi- 
component KP hierarchy [96] and therefore by the following Baker-Akhiezer function 
system: 



5str^*(t;C) = Q(t;5)^(t;C)- 



(2.1) 
(2.2) 



Here the operator P{t;d) and Q{t;d) are p-th and g-th order differential operators in 
d = gstrdt, respectively, which satisfy the Douglas (string) equation \21\ : 



[P{t;d),Q{t;d)] =g,trlk 



(2.3) 



Critical points in the multi-cut two-matrix models are characterized by these Lax oper- 
ators and explicitly obtained in [16] with their critical potentials. There are two kinds 
of interesting critical points: the Z^-symmetric critical points and fractional-superstring 
critical points. A brief summary of the corresponding Baker-Akhiezer function system is 
following]^ 

1. The Zfc-symmetric critical points are characterized by the following k x k Lax 
operators [46j: 

P{t;d) = Td^ + J2uL^'^\t)d'', Q{t-d) = T-^d^ + ^Ui^^^\t)d'', (2.4) 

n=0 



n=0 

with the shift matrix F 



/O 1 

1 



V 



1 
0/ 



(2.5) 



and the k x k matrix- valued real coefficients U^''^\t) and Un^''^\t) which satisfy 

/o * \ /o A 

* * 

[/f^^)(t)= •.. •.. , f/f'=«)(i)= , (2.6) 



* 

0/ 



V 



* 

* oy 



'For the derivation of these systems from the multi-cut two-matrix models, see 
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as a result of the symmetry of the critical points. Macroscopic loop amplitudes 
(i.e. off-critical resolvent amplitudes with t 7^ 0) in this kind of critical points are 
also obtained in |16] with the Daul-Kazakov-Kostov prescription |30] and expressed 
as the Jacobi polynomials or the third and fourth Chebyshev polynomials. In 
particular, the amplitudes in the the fc-cut (1, 1) critical points are given as the 
eigenvalues of the Lax operators Eq. (12. 4p in the weak coupling limit Qstv — > 00 



P(t;a)^diag^ti(piLeai(^;^; 
Q(t;9)^diag?=ifQSssicai 



diagj^i 



u^-^ x(z] 



[t-z] 



diag' 



with 



=t\l [z-c)\z-hY \ y{z) = t\l {z-cY \z - h)^ 
and Q = cl + h {k — I) and the dimensionless variable z = Qstrt'^dt. 



(2.7) 



(2i 



2. The fractional-superstring critical points |13] are characterized by the following two 
kinds of Lax operators [16]: The first kind is given as 



p-1 



q-l 



P(t;9) =r9P + ^t/f^^)(t)a", Q(t;a) = r9^ + ^f/f*«)(t)5". (2.9) 



n=0 



n=0 



These Lax operators are derived from the cj^/^-rotated critical potentials. The 
second kind is given as 



p-1 



Pit- d) = r^"^"^') d^ + Y^ ui^^^\t) a", 



n=0 



(2.10) 



n=0 



with the matrix r*^''*"^'^ 



/O 



■p(real) 



1 

1 



1 



(2.11) 



These Lax operators are derived from the real critical potentials. In both cases, 
all the k X k matrix- valued coefficients ulf''^\t) and uf^'^\t) (and Un^''^\t) and 
are real functions. The macroscopic loop amplitudes in each case are 
obtained and given by the deformed Chebyshev functions [43|. 

In this paper, for the sake of simplicity, we concentrate on the p = 1 cases of the 
Zfc-symmetric critical points. With this choice of critical points, the Lax operator P{t; d) 
becomes 



Pit; d) = Td + Hit), 



(2.12) 
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In this paper, the equahty ~ means that they are equal up to some similarity transformation. 
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and the Baker- Akhiezer function for the eigenvalue problem of the operator P{t;d), 
Eq. (12. ip . is rewritten as 

gstr^^'^it; = vit- ^{t; = r-i [c - H{t)] ^{t; c), (2.13) 

and therefore Eq. (12 ■2p is also rewritten as a k x k matrix polynomial operator in (: 
9st.^^^ = Q{t; C) ^(t; C) = Qit; d) ^{t; C), Q(t; C) = E Q-n{t)C-'- (2.14) 

n=l 

Here we define r as 

r = g + l>0, (2.15) 

which is referred to as the Poincare index in literature. The advantage of this formulation 
is that the pair of Lax operators {P{t;d),Q{t;d)) becomes a pair of the polynomial 
operators (V (t-X) , Q{t-X)) , and the system can be expressed as an /c x first order 
ordinary differential equation (ODE) system. These systems are called the Zakharov- 
Shabat eigenvalue problem [97] or AKNS hierarchy [98] in literature. Note that the 
Douglas equation becomes 

[P{t;d),Q{t;d)] =g,,Jk ^ [9stA ' Qit; 0, 9stA - V{t; ()] =0, (2.16) 

in terms of these Lax operators. 

This ODE system Eq. (I2.14p has the k independent order k column vector solutions 
"^^^Kt; C), (j = 1, 2, ■ ■ ■ , k), and we here use the following matrix solution notation: 

^{t-0^[^^'\t;0,... ,¥'^\t-X)y (2.17) 

As in the usual ODE, we consider formal expansion around C — t- oo. However the 
point C ^ oo is an irregular singularity and the formal series expansion around this 
irregular point in general does not converge absolutely. Up to proper redefinition of the 
k independent solutions, the formal series expansion of the solutions around — oo is 
given as 

W«)->'(«;0ei*«-[4 + f;^]xexp[-L(^„luC- E ^)'. 

n=l ^ ^^^^ m=-r,^0 ^ 

(2.18) 

The coefficient matrices are obtained from the recursive equations, 

n+r 

= -ngstrYnit) + ^ ^m(t) iPn-mit) - Qn-rn{t) Fm(t) , (?2 = -r, -r + 1 , ■ • ■ ) . 
m=0 

(2.19) 

For convenience, we extend the indices of the coefficient matrices: 

Fo(t) = 4, Y^{t) = {n<0), ^^{t) = QUt) = (m < -r), (2.20) 
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' ^^^^ m=-r,^0 ^ 



and impose the following constraints on Yn{t) and (pn{t)- 

k 

[r', Mt)] =0, Yl i^M = 0' (/ = 0, 1, ■■■,/.- 1). (2.21) 

i=l 

This recursive relation then can be solved uniquely and all the expansion coefficient are 
written with the coefficient matrix- valued function H{t) in Eq. f l2.12p . 

On the other hand, it is also convenient to use a diagonal basis, "^asymit] C); which is 
defined by 

It c 

= U^-^asym{t;C)U, (2.22) 

where the matrix U is given as 

Uji = TU = Un, (2.23) 

with = diag(l, w^, ■ ■ ■ , w'^~^) and co = e^'^*/^. Since this is a similarity transforma- 
tion, the coefficients also satisfy the same recursive relation f l2.19p . In this basis, the 
function ip{t; () is a diagonalized matrix and we write its eigenvalues as 

^{t; = diag(^«(t; C), ■ " " , ^^'\t; C))- (2-24) 

The vector components of the formal series, '^asym = {^asym, ■ ■ ■ , '^itsym) , is given as 

^l^lmit; C) = Y('\t; C) e^^'''^'-'^\ (j = 1, 2, ■ ■ ■ , fc), (2.25) 

with Y{t;C) = - 

Although the above formal solutions are formal series around the irregular singularity, 
they are related to the exact analytic solutions of the ODE system, \l/(t; (), in the sense 
of asymptotic expansion: 

^{t;C) ^ ^asym{t;OC, (2.26) 

asym 

in some specific angular domain [99] : 

C -> oo G D{a,b) = {C G C; a < arg(C) < b}. (2.27) 

An example of the angular domain is shown in Fig. [2]-a. Here C is a proper coefficient 
matrix, and the meaning of asymptotic expansion is following: 

Definition 1 (asymptotic expansion) For a holomorphic function /(C) ? an asymp- 
totic expansion of f{() in a domain D{a,b) is defined as a formal series ^„/nC~" such 
that there exists a constant -B^\ G M which satisfies 

/(O - E ^1 < ^ ^ ^(«'^) n {C e C; Id > R] (2.28) 



n=—r 



for each integer N = — r, — r + 1, ■ ■ ■ and sufficiently large i? G M. This is written as 

oo „ 

/(O ^ EtS' C^ooeD{a,b). (2.29) 

n=—r 

The maximal angular domains are called Stokes sectors. 
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2.2 General facts on Stokes phenomena in the ODE system 

In this subsection, in order to understand the asymptotic expansion Eqs. fl2.22p and 
( ]2.26p . we review some general theorem about the asympototic expansions and Stokes 
phenomena in the general k x k ODE systems, 

5str^5 (t; C) = [Q-r C'-' + Q-r+lit) C"' + " " " Q-i{t)] nt; C) 

^Q(t;C)$(t;C). (2.30) 

Note that proof of the theorems appearing in this subsection can be found in ^1] and 
references therein. For sake of simplicity, we assume 

= diag(Ai,A2,--- Ai-Aj^O, A 7^ 0, (i, j = 1, 2, ■ ■ ■ , fc). 

(2.31) 

Therefore, the exponents Eq. fl2.24p are expressed as 

^(t;C) = ^o(t)lnC- E ^ = -Q-rC + --- , (2.32) 

n=—r,nj^O 

and v?^^ = Ai {i = 1,2, ■■■ , k) also satisfies fl2.3ip . 

The meaning of the asymptotic expansion Eq. f l2.22p is that basically we ignore rela- 
tively small exponents. One takes some (small enough) anglular domain D{a,e^'^a) then 
compares the relative magnitudes around ( ^ 00, for example, 

|e^«^'(*;C)| < |e^*^'^'(*^^)| <...< |e^*^'^-'(*^^)|, ( ^ 00 e D{a,e'%). (2.33) 

Then one can obtain the following equality under the asymptotic expansion: 

e^'^'^nt;C) + ^^e^'^i'(*;C) ^ e^«2){t;C)^ C ^ 00 G Z^(a, e^^a). (2.34) 

asym 

That is, the smaller exponents become practically invisible in view of the asymptotic 
expansion. Our interest is to identify angles of ( where the the exponents, exp((y9*^-'^(C)) 
{i = 1,2, ■ ■ ■ , k), change the relative magnitudes around C — )■ 00. This leads to the 
concept of Stokes lines: 

Definition 2 (Stokes lines) With the assumption ^2. 31\) . Stokes lines SL^y in this ODE 
system are defined for each pair of (j, /) as 

2r-l 

SL,, , ^ {C G C; Re[(^^] - (^«)Cl = o} = IJ SlJJ, (2.35) 

n=0 

which consists of 2r semi-infinite lines, SL^"'' {n = 0,1, ■ ■ ■ , 2r — 1) . The set of lines, SL, 
denotes a set of whole Stokes lines, SL = (J^ ^SL^y. 

An example of Stokes lines SLj^; is shown in Fig. |2]-b. In particular, if the angular 
domain D{a,b) of the asymptotic expansion includes a Stokes line, one cannot neglect 
the exponents as it happens in Eq. fl2.34p . This leads to the following definition of Stokes 
sectors: 
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Definition 3 (Stokes sectors) A Stokes sector D in the ODE system is an angular 
domain, D = D{a,b), with angles a and b such that for each pair of there exist a 
unique Stokes line SL^"^'''' which satisfies, 

^hff^ <Z D = D{a,h), (2.36) 

that is, except for this line SL^"-''''' there is no other line SL^ /'''' SL^"^''''j which runs 
inside the domain, D. 

An example of the Stokes sectors (the 3-cut (1, 1) critical point) is shown in Fig. |2]-b. 




Figure 2: a) An angular domain of I?(a, b). b) Stokes lines and Stokes sectors. This is the 3-cut (1, 1) 
critical points. An example of Stokes sectors is also shown. In this critical point, there are three kinds 
of the Stokes lines SL^j-, (i, j) — (1, 2), (2, 3), (3, 1). Stokes sectors includes one and only one Stokes line 
of each kind. 

Actually the definition of the Stokes sectors results in the following theorem [99] : 

Theorem 1 For a given Stokes sector D, any solutions to the ODE system \l/(t; Q) has 
the following asymptotic expansion: 

^{t-Q ^ ^asym{t-X)C, C OO G A (2.37) 

asym 

with a matrix C. Furthermore, the coefficient matrix C (i.e. asymptotic expansion) is 
unique in the Stokes sector D. 

This uniqueness enables us to define the following unique solution in a Stokes sector D: 

Definition 4 (Canonical solution) If the solution to the ODE system, can{t'-, C) j has 
the asymptotic expansion with C = Ik in a Stokes sector D, 

$ea„(t;C) ^ ^asymit-X), C ^ OO G (2.38) 

asym 

this solution is called the canonical solution in the Stokes sector D. 



11 



This theorem on the other hand means that the asymptotic expansion is not unique 
if one chooses some angular domain D' narrower than Stokes sectors. In particular, as is 
shown in Fig. |3l the intersection of two different Stokes sectors Di and D2 is generally 
narrower than Stokes sectors, and therefore there appears difference between the canonical 
solutions \E'j(t; C) of each sector Di{i = 1, 2): 

$2(t;C) = ^i(t;C)5, AnD2 7^0. (2.39) 

This k X k matrix S which expresses the difference between \E'i(t; C) and '^2{t] C) is called 
a Stokes matrix in the intersection Di (1 D2. This indicates that solutions in the ODE 
system generally have different asymptotic expansion in different Stokes sectors. This 
analytic behavior of the solutions is referred to as the Stokes phenomenon in the ODE 
system. 

^2{t;0 ^ $asym(t;C) 

(C ^ 00 G D2) 

C) = ^l(t;C)5 ^ ^asym(t;C) 
asym 

(c ^ 00 G Di n D2) 

C) ^ ^asym(t;C) 
asym 

(C ^ 00 G Di) 




Figure 3: Explanation of Stokes phenomenon in ODE systems. For given two Stokes sectors, their 
canonical solutions are generally different by a Stokes matrix, S in the intersection Di O D2- This 
behavior of analytic functions is called Stokes phenomenon. 



A direct calculation shows that the Stokes matrices do not depend on (, and further- 
more, they do not depend on the deformation parameter t either (as in (I2.13P ): 

dS dS / , X 

^ = ^ = 0. 2.40) 

dC dt ^ ^ 

This means that the Stokes matrices are understood as integration constants for the 
evolution system in the t space. Therefore, these integrable deformations in the original 
multi-component KP hierarchy are also called isomonodromy deformation system [84| . 
This also leads us to the concept of inverse monodromy approach [8ll[85], which is also 
briefly reviewed in Section [51 

Components of Stokes matrices satisfy the following theorem (See [91], for example): 



Theorem 2 (Stokes multipliers) For given two Stokes sectors Di and D2 (Di r\D2 ^ 
components of their Stokes matrices, i.e. Stokes multipliers, S = (sj.j), satisfy 

s,,, = l {j = l,2,--- ,k), (2.41) 
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and Si J {i 7^ j) can take non-zero values only when the exponents satisfy 

Re[ip^llC] < Re[ip^ilC] for all angular range of C ^ 00 E Di n D2 ^ (2.42) 
In particular, these Stokes multipliers, Sij {i 7^ j), are called "non-trivial" . 
In this paper, we often refer to these facts about Stokes phenomena in ODE systems. 

2.3 Stokes phenomena in the two-cut case 

In this subsection, we speciahze the general consideration to the two-cut (1,2) critical 
point. 



2.3.1 The ODE system and asymptotic expansions in the two-cut case 
The string equation in this system is known as the Painleve II equation [351136] . 



9str 



/-f + 2t/ = 0. 



(2.43) 



which is equivalent to the following ODE system in ( (Eq. fl2.45p ) with its isomonodromy 
deformations in t (Eq. (I236D)0 



9str-^^{t; C) = [c^aC - ^1 fit)] ^{t; 0- 



1 \ i 

-f^ + lijaz- g^tr^a^f ^(t;C), 



(2.45) 
(2.46) 



Since this 2x2 first-order ODE system has two independent column vector solutions 
\E'^^^(t; Q and '^'^'^\t] we use the matrix notation for the solutions: 



^>{t-0={^^'\t-QM'\t;0)- 



(2.47) 



At the point C — > cxd, the ODE has an irregular singularity (of the Poincare order 3) 
and the formal expansion of the solutions fl2.22p is given as 



asymiC;t)= [h + ^a2f (t) + 0(1 / 



exp 



— (\cr,C' + IiasC + 0(1/0 

5'str ^3 



= y(t;C)e 



(2.48) 



This can be obtained with the recursion relation Eq. (12.1 9p (see also in Appendix IB.2p . 
Note that the exponent ip{t; () is a diagonal matrix which satisfies (f{t; () oc cxs, and then 
each vector solution ^aiym(t; C) (^ = 1)2) has different exponents: 



with 



r(t; C) = (r«(t; C),Y'^'\t; C)), ^(t; C) = diag(^«(t; C), ^^'\t; ()). 



(2.49) 



(2.50) 



^-'^In the later discussion (from Section[3]), we also define a different basis: ^(t; C) = U'i>{t; C)t^^ with 

[/cTgf/t = fji, [/cTiC/t = -CT3, Ua2U'' = CT2. (2.44) 

This basis naturally appears in the matrix-model calculations and is more suitable to read the Hermiticity 
of the multi-cut matrix models 1461. 
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2.3.2 Stokes sectors and Stokes matrices 

In this case, there is only one kind of the Stokes lines SLi^2 which is given by fl2.35p as 
Therefore, Stokes sectors Dn are given as 



(n = 0,l,--- ,5). 



(2.51^ 



(n = 0,l,--- ,5), Do = D(-^,^). 



0, 



(2.52) 



This is shown in Fig. HI The canonical solution on the Stokes sector Dn is denoted by 
\E'„(t; Q. The Stokes matrices Sn are now defined as 



Sn = ^-\t-C)^n+i{t-X), (n = 0,l,--- ,5), 
and therefore components of the Stokes matrices are read as 



D2n n D2n+l 



s. 



2n 



1 

S2n 1 



D2n+l n D2n+2 '■ S2n+l — ( g ^^^^ ; ; 



C ^ oo , 



(2.53) 



3^W(*;C)| < |e^<^Ht;C)| C^oo). (2.54) 






(b) 



(c) 



Figure 4: a) Stokes lines in the two-cut (1,2) case, b) Stokes sectors of D^.D^ and -D4. c) Stokes 
sectors of Di,Dt, and D^. 



2.3.3 Three basic constraints on the Stokes multipliers 

The Stokes multipliers satisfy three constraints from the symmetry of the original ODE 
system. 



Z2 symmetry constraint This symmetry originates from the Z2 symmetry of the 
matrix model. That is, this is the reflection symmetry M — )■ —M of the one-matrix 
models: 



^^g-7Vtry(M)^ V^(-M) = V{M). 



(2.55) 
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In terms of the ODE system, this symmetry is expressed by the reflection of C — ^ — C- 

9str^^^^^ = [-Qit;-C)] ^(t;-C) = [^iQ(t;C)^i] 5(t;-C), 

9str^^^^ = [nt; -0] ^(t; -C) = hnt; C)ai] ^(t; -()• (2.56) 
Therefore, each canonical solution is mapped to another canonical solution as: 

-0^1 = $„+3(t; C), (n = 0, 1, ■ ■ ■ , 5), (2.57) 
and the Stokes matrices are mapped as 

Sn+3 = CTlSnC^l, 3^+3 = Sn, (n = 0, 1, ■ ■ ■ , 5) . (2.58) 

Consequently there are only three independent Stokes multipliers. 

So = S3 = a, si = S4 = /3, S2 = S5 = 7. (2.59) 

Hermiticity constraint This originates from Hermiticity of the matrix models. In 
the two-cut cases, they are studied in [39l[69]. This symmetry is expressed a^ 

Q*it; C) = Qit; C), V*it; C) = Vit; C). (2.60) 

Therefore, each canonical solution is mapped to another canonical solution as: 

%{t; C) = ^7-n{t; C), (n = 0, 1, ■ ■ ■ , 5), (2.61) 

and the Stokes matrices are mapped as 

S: = S^\, s: + S6-n = 0, (n = 0,l,--- ,5). (2.62) 

This reduces three independent Stokes multipliers a, /3 and 7 to be two real parameters: 

a* + a = 0, /3* + 7 = 0. (2.63) 

Monodromy free constraint The last constraint is the requirement that the solutions 
to the ODE system are single-valued functions. Note that the presence of non-trivial 
monodromy in the context of matrix models corresponds to introducing background RR 
flux and/or DO-branes in OA string background. That is, the system becomes like the 
complex matrix models P2 | |88 |[T0T] . This constraint for the single-valued solutions is 
expressed as 

$„(t; C) = $n(t; e'"C) = 5n+6(t; C), (2.64) 

therefore 

^0^1^2^3^4^5 = /2, (2.65) 

^^Note that we use the following convention of complex conjugation in this paper: [/(C)]* = /*(C*) — 
E„ /*C*, with a function /(C) = E„ fnC- 
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which results in 

S0 + S1 + S2 + S0S1S3 = a(l - 1/31') +/3-/3* = 0. (2.66) 

Taking all constraints Eqs. fl2.59p . (I2.63p . fl2.66p into consideration, we find that the 
Stokes multipliers have two real degrees of freedom, say /3. Since the Painleve equation 
II equation, Eq. (12.43^ . is the second order ODE system. These two parameters are the 
non-perturbative ambiguity of the system. 

As is mentioned in Introduction, among these Stokes multipliers satisfying the alge- 
braic relation f l2.66p . there is a special value which realizes the perturbative behavior (in 
t ±00) of the matrix models argued from the physical point of views J42j. This special 
value is given as 

a = 0, /3 = ±1, (2.67) 

and corresponds to the Hastings-McLeod solution in the Painleve II equation [89]. From 
the mathematical point of view, this solution also has a good analytic behavior along 
the real isomonodromy parameter (cosmological constant) t [89 |ll00j . From this two-cut 
example, we generally expect that there is a special class of solutions of the Stokes multi- 
pliers which corresponds to the physical D-instantion chemical potentials. To generalize 
the solutions to the cases of arbitrary number of cuts, it is natural to ask the following 
question: what is the physical requirements which specify the above multipliers? This is 
also related to the issue cited by [6ll|63]: What is the boundary condition in continuum 
formulations which can fix the D-instanton chemical potentials in the matrix models? 
Our procedure (discussed in Section H] and Section [5]) gives an answer to the question. In 
Section 14.2.11 and then in Section ISTTl we will see that our physical requirements correctly 
choose this particular parametrization Eq. fl2.67p of the Stokes multipliers. 



3 Stokes phenomena in the multi-cut cases 

In this section, we develop general framework for Stokes phenomena in the general multi- 
cut critical points, and show explicitly how the actual systems can be controlled. Key 
information is provided by profile of dominant exponents (Theorem [3]), and with this 
terminology we propose a systematic way to read the non-trivial Stokes multipliers (The- 
orem H]). Since the following discussions are valid in general k x k ODE systems of the 
following type: 

^[51|0 = (r-7^-i + ...)v[;(t;0, g.c.d.(A;,7) = l, (3.1) 

we here develop the general framework without restricting to the symmetry (7 = r). 
The restriction to the symmetric cases only appear in Section 13.31 

3.1 Stokes lines and Stokes sectors 

First we focus on the Stokes lines, 

SL,-,: Re[(^^]-^('i)Cl =0, (3.2) 
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and the resulting Stokes sectors fj2.36l) . The leading coefficient of the exponents, ^p^ll, 
which we consider here is given a; 
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<^(f) = w-tO'-D. (3.3) 

Consequently, the conditions on the Stokes lines (in terms of angle, ( = |C|c*^) are 
expressed as 

Re[{^^l - y^^De^'-^] = 2sin(r^^ - 7r ^^^'+/~^^ ) sin(7r^^^:^) . (3.4) 

First of all, if there is a pair of (j, I) such that 

7(J - e kZ, (3.5) 

then the condition fl2.3ip does not satisfy. This means that the highest exponents degen- 
erate (v?^] - ip^ll) C = 0. In this case, we consider the next leading Stokes lines, 

M{^-li-y'^-l+i)C-']=0, (3.6) 
or more generally we consider the following Stokes lines0 

Definition 5 (General Stokes lines) The general Stokes lines GSLj^; in this ODE sys- 
tem are defined for each pair of (j, /) as 

2r-l 

GSL,-z ^ {C G C; Re[^(^)(t; C) - ^^'\t; ()] = o} = |J GSL^J^^ , (3.7) 

n=0 

which consists of 2r semi-infinite lines, GSL^"^ (n = 0, 1, ■ ■ ■ , 2r — 1). The set of lines, 
GSL, denotes a set of whole (general) Stokes lines, GSL = IJ^.^GSLj^;. 

The situations (13. 5p are also interesting critical points in the multi-cut matrix models, 
however here for sake of simplicity, we concentrate on the following cases, 

g.c.d.(A;,7) =1, (3.8) 

because Eq. (13. 5p becomes trivial in this case: 

7(j -I) ekZ ^ j -I ekZ. (3.9) 

Therefore Eq. (13.41) gives the angle O^"'^ for the Stokes hues SL^^; as 

a ziW fcn + 7(j + /-2) 

^ = = ] n e Z. (3.10) 

From this formula, one can read several basic information about the Stokes lines. An 
example of Stokes lines (3-cut (1, 1) case) is shown in Fig. [2]-b. For later convenience, we 
introduce the following terminology: 



""^^Note that the cases of our interest in the later sections are the Zfc-symmetric critical points, and as 
one can see in Appendix [Bl the cases are given by 7 = r. Also for future reference, we note that the 
fractional-superstring cases are given by 7 = r — 2. 

^"^Thc physical interpretation of these general Stokes lines is the positions of eigenvalues in the matrix 
models. This viewpoint is also essential in this paper and discussed in Section [ 
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Definition 6 (Segments) Angular domains in between two Stokes lines which do not 
include any Stokes lines are called segments. 

In our present cases with a coprime (fc,7), there are 2rk distinct segments 5Dn {n = 
0, 1, ■ ■ ■ , 2rk — 1) given as 

6Dn = D{n6e -6e,n6e), (n = 0,1,- ■ ■ ,2rk - 1; 66 = (3.11) 

which can fill the complex plane C, 

2rfc-l 

IJ Il\ = C, SDm n 6Dm' = {m^m'). (3.12) 

n=0 

According to the definition of Stokes sectors, Eq. f l2.36p . we define the following most 
basic Stokes sectors, 

Definition 7 (Fine Stokes sectors/matrices) The following angular domains Dn 

D„ = e™^^Do, Do = D{-6e,k6e), (n = 0, 1, ■ ■ ■ , 2rA; - 1), (3.13) 

are Stokes sectors of a coprime {k,r) system with k > 3, which are referred to as fine 
Stokes sectors. The canonical solution of the fine Stokes sector Dn is denoted as \l/„(t; C) 
and the corresponding Stokes matrices Sn are given as 

^n+i{t;C) = ^n{t;C)Sn, (3.14) 
which is referred to as (fine) Stokes matrices. 

Here we also define the other two kinds of Stokes sectors/matrices: First we define Stokes 
sectors/matrices which respect to the symmetry of the multi-cut matrix models: 

Definition 8 (Symmetric Stokes sectors/matrices) The following subset of the fine 
Stokes sectors, 

D^nr. (n = 0,l,--- ,A;-1), (3.15) 
are referred to as symmetric Stokes sectors^^ and the corresponding Stokes matrices 

ri(sym) 

i(sym) 



S2rn = ^2rn(^i C) 2r{n+l){t', C) = S2rn ' 5*2™+! " " " S2r{n+l)-l- (3.16) 

are referred to as symmetric Stokes matrices. 

Next we define the following economical Stokes sectors/matrices: 



^^Note that this definition is not enough for the k — 3,r = 2 case. In these cases, we employ a modified 
version of the Stokes sectors, for example, Dnr- 
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Definition 9 (Coarse Stokes sectors/matrices) The following subset of the fine Stokes 
sectors, 

Dr^k, (n = 0,l,--- ,2r-l), (3.17) 

are referred to as coarse Stokes sectors, and the corresponding Stokes matrices S^^jj. are 
written as 

^(coa.) ^ Q \E'(„+i)fc(t; C) = Snk ■ Snk+l " " " (3.18) 

are referred to as coarse Stokes matrices. 

Coarse Stokes sectors are most often used in the literature. However, in the follow- 
ing discussions, one will see that the fine Stokes matrices are more convenient for our 
calculations. 



3.2 Stokes multipliers from the profile of dominant exponents 

In principle, one can use Theorem[2]to read the non-trivial (or non-zero) Stokes multipliers 
in each specific case. That is the problem of finding which components can take non- 
zero value in Stokes matrices. However, practically in general, it is tedious to use this 
standard way to read the non-trivial multipliers, especially in the higher k x k system 
with higher Poincare index r. The purpose of this section is therefore to point out an 
interesting connection between the non-trivial Stokes multipliers and profile of dominant 
exponents which we develop in this subsection (Theorem 12] and H]). An important thing 
in this procedure is that these results make it easy to put data of the Stokes multipliers 
in computer, for example, in Mathematica program. 

Since there is no Stokes line in the segments defined in Eq. ( 13. lip , one can define the 
following ordered set Ji of indices ji^i: 

Ji = [ Ji,i\ji,2 I ■■■ \3i,k ] GN^ (3.19) 
which describes the profile of dominant exponents in the segment Di 

Re[^^';^'>C] < Re[^^^;'\^] <■■■< Refy.^^-'^^Cl , C e SD^. (3.20) 

This sequence of numbers, J' = { J^jflg, is referred to as profile of dominant exponents. 
Here we express the profile J7 as follows: 



j2rA;-l,l 


j2rk-l,2 




j2rk-l,k 










Jl,l 


jl,2 




jl,k 


J 0,1 


jo,2 




jo,k 



Note that the ordering of indices in the vertical direction is different from the usual 
matrix, and that elements are periodic in the index I, Ji = Ji+2rk- An example (3-cut 
(1, 1) critical point) and the relation to the ( plane are shown in Fig. |5l 

The non-trivial problem for the profile is then how to fill the numbers in the profiles. 
We found the following simple answer: 
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(a) (b) 

Figure 5: The two expressions for the profile of dominant exponents in the 3-cut (1,1) critical point. 

Re[i^^2'^] < Re[</5^2^^] < ^'AH^^--f'\ is expressed as < j/_2 < j;,3- a) The profile in the C plane, b) 
The profile J in the table. In the same way, the dominance is expressed as [ji,i|i/,2b/,3] 



Theorem 3 (General cornponents) The general components ji^n of the profile J with 
g.c.d. (k,^) = 1 are given atJ 



+ 



\k-\-l+n 



k — n + 1 



mi, 



mod /c, 



(3.22) 



where mi is obtained by the Euclidean algorithm of kni + 7ml = 1. 

Some comments on this formula are in order: 

• In any segment profile J; = [ ji^i | ji^2 | ' ' ' | ji,k ], a pair of indices {i,j) which 
change their relative dominance at angle 6 = 166 are always next to each other, and 
they satisfy the following sum rule0 



i+j 



nil 



Imi] 



mod k. 



(3.24) 



with an integer mi which is obtained by the Euclidean algorithm of kni + '-ymi = 
I. Therefore in particular, we represent these pairs as {i\j) in the profile (See 
Eqs. (ESSD). 



Theorem E] can be recursively shown by using the sum rules Eq. fl3.24p . and its 
initial conditions jo,fc = ji,fc = 1 = jk,i = jfe+1,1 which can be easily checked. 



""^^In this paper, we use the floor-function notation for the gauss symbol, [aj , which means the largest 
integer less than or equal to a. 



^^Note that the condition for ~ ISO is given as 



fcn + 7(1 + 7 - 2) , , „x , 

^\ '-TT ^ —TT ^ kn + j{i+j-2) = l. 

rk rk 



(3.23) 



This means that, for a given Z, find a pair such that there exists an integer n. 
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The trajectories of the indices, for instance a and 6, are given as follows: 
is odd A; is even 



J 



















a 






b 






a 










b 




a 












b 




a 
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b 



















J 
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a 












b 


a 












b 




a 








b 








a 




b 





















(3.25) 



The above formula is not applicable in the case of g.c.d (/c,7) 7^ 1, since some 
exponents degenerate, but one example of this kind is shown in Appendix [Fl 



Here also two examples of the profiles J7fc,r are shown for the case of {k, r) 
(5,2) with 7 = r (Z^ symmetry condition): 
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(3.26) 



^ Jo 

One can observe that there is a 2k periodicity, ji+2k,n = ji,n, or more precisely, a reflection 
by step k, ji^n = ji+k,k-n+i- 

Next we demonstrate how to read the non-trivial Stokes multipliers in some examples, 
and see the general rule. In the case of (r, k; 7) = (2, 5; 2) and its symmetric Stokes matrix 
Sq^'^\ one first sees the dominance profile in the domain DqC] D4, 



^ J5 
J4 



(3.27) 
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and reads the ordering of magnitude: 



(2)>(5),(4),(3),(1), (5)>(3),(1), 
This results in the following symmetric Stokes multipliers: 



(4)>(3),(1). 



(3.28) 
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(3.29) 



In the same way, for the calculation of the fine Stokes matrix 5*0, one first sees the 
dominance profile in the domain Dq fl Di, 
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(3.30) 



and reads the ordering of magnitude: 

(4) > (3) 

This results in the Stokes multipliers: 

/I 



(2) > (5). 



(3.31) 
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(3.32) 



These are the standard way of reading the multipliers! 

However, one may notice that there is a relation between indices of non-zero Stokes 
multipliers So,j,j in the Stokes matrix Sq and the dominance-changing pairs in the 
profile Jq: 



S0,3,4, So,5,2 ^ (2|5), (4|3) G Jo 

We claim that this observation is generally true: 



[ (2 I 5) I (4 I 3) I 1 ] 



(3.33) 



Theorem 4 (Stokes multipliers from the profiles) The non-zero Stokes multipliers 
in the fine Stokes matrix Si have a correspondence with dominance- changing pairs {j\i) 
in the profile Ji as follows: 



{i 7^ j) can take non-zero value ^ {j\i) G J/. 
Note that the orderings of indices {i\j) and sij^i are opposite i ^ j. 



(3.34) 



-f(coa) 



18 Pj-oitl this procedure, one may notice that the simplest choice is the coarse Stokes sectors S]^^' , 
because intersections have the definite order of magnitude: DoflDsD [l|3|4|5|2] and the 
number of Stokes matrices is the smahest. This is the main reason why the coarse Stokes sectors are 
often used in hterature. However, we wih see that the coarse Stokes muhiphers are not suitable for 
general formula of higher k and r at least in the symmetric critical points. 
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A proof is easy if one notices that intersections of fine Stokes sectors D„ fl -D^+i are 
always a half of the period of Stokes line formula Eq. (13.41) . The other Stokes matrices, 
say S'n and Sn°^^ , are written as a product of the fine Stokes matrices Sn (as in (I3.16P 
and (I3.18P ). For instance, 
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As one can see from these special examples, the Stokes multipliers are always related as 
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(3.35) 
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and one can then show that the number of independent Stokes multipliers in each Stokes 
matrix Sq'^^'^ and 3^^°^^ is the same and is supplied by the fine Stokes matrices. Conse- 
quently, the same statement also holds for these different kinds of Stokes multipliers: For 
example. 



4r(™\) ij 7^ i) '^^^ ^^^^ non-zero value 

^ ij\i)eJn, n = 2r(/-l),2r(Z-l) + l,--- ,2H-1. 



(3.37) 



3.3 Three basic constraints on the Stokes matrices 

Finally we show the three basic constraints on the Stokes multipliers, which result from 
the detail analysis of (the Z^-symmetric) critical points in the multi-cut two-matrix mod- 
els [IS] and also which provide a natural extension of the two-cut cases (See Section 1^75]) . 
We should note that the conditions from the matrix models are given in the F-basis (or 
the matrix-model basis) \E'(t;C) and the Stokes matrices are defined in the r2-basis (the 
diagonal basis) \l/(t; (), and they are related by a unitary transformation (See Eq. (12.221) 
and Eq. fl2:23|) ). 

Zfc symmetry condition This condition is from the symmetry in the multi-cut 
two-matrix models [46j and generally expressed a^ 

uj~^Q{t;uj-\) =n-^ Q{t;C)n, V{t;io-\) = n-^V{t;C)n, (3.38) 
^^This is a direct consequence of Eqs. (I2.6p . 
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with Q ^E'j = uEi^i^i. The constraint on the Stokes matrices are then obtained as 

Sn+2r = T-'SnT, {n = 0, 1, ■ ■ ■ , 2rk - 1) (3.39) 
for the fine Stokes matrices S'„. A note for the derivation is following: 



1. Because of the condition Eq. f l3.38p . the canonical solution \E'„(t;C) for a Stokes 
sector Dn satisfies 



(3.40) 



and therefore one obtains 

^„+2.(t; C) = [^^*n(t; UJ~^Q^-^\ ^asym(t; C), C OO E Dn+2r = D^. 

(3.41) 



asym 



2. By translating this relation into the Q-basis (diagonal basis), 

^n+2r(t; C) = fi^n(t; C^^'C)^"' ^ ^n+2r(t; C) = T-'^nit, UJ'\)T , (3.42) 

with U^^QU = F"^, one obtains the relation of the Stokes matrices: 

Sn+2r = KUrit; C)^n+2.+l(t; C) = ^-'K^t; C)^n+l(t; or = T-'SnT. (3.43) 

This condition means that only the first 2r Stokes matrices Sn{n = 0,1, ■■ ■ , 2r — 1) are 
independent. Therefore, we use the first 2r dominance profiles to identify the non-trivial 
Stokes multipliers: 
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(3.44) 



Here we show an examples ofA; = 5,r = 7 = 2: 
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(3.45) 



Hermiticity condition This condition is from the hermiticity of the multi-cut matrix 
models ^46j and generally expressed as 



Q*{t;C) = Q{t;C)- 

The constraints on the Stokes matrices are then obtained as 



(n = 0, 1,-- - ,2rk-l) 



(3.46) 



(3.47) 



for the fine Stokes matrices Sn- A note for the derivation is following: 
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1. Because of the condition Eq. fl3.46p . the canonical solution \E'„()f:;^) for a Stokes 
sector Dn satisfies 

g.,r^f^ = Q{t-X)Kit-X), (3.48) 

and therefore one obtains 

^'(2.-l)fc+l-„(t;C) = ^:(t;C) ^ ^asym(t;C), C ^ /^(2r-l)fc+l-n = (3.49) 

asym 

2. By translating this relation into the Jl-basis (diagonal basis), 

Vl/:(t; C) = ^i2r-l)k+l-n{t; C) ^ = ^(2.-l);t+l-n(t; C) 

= Ar $(2._i)fc+i_„(t; C) r-^A, (3.50) 
with U* = and U"^ = AT, one obtains the relation of the Stokes matrices 

s: = [^-\t-o^n+i{t;OY = Ar[$^2l-i)m-n(^;C)$(2.-i).-n(t;C)]r-^A 

= Ar [^i,U),_^it; C)$(2r-i).+i-n(t; 0] A 
= Ar^(-^_y,_r-iA. (3.51) 

Monodromy free condition If the formal expansion satisfies (/Jq = (discussed in 
Appendix [B]) , then the canonical solutions are the single- valued functions: 

^n{t; C) = ^n{t; e'^X) = ^n+2kr{t; C), (3.52) 
therefore the Stokes matrices satisfy 

ri ri ri ri(coa) (~i(coa) (~i(coa) 

Do ■ Oi ■ ■ ■ D2rfc-1 — Oq ■ D;^ ■ ■ ■ 0^(2r-l) 

_ c(sym) r.(sym) Q(sym) _ j rn\ 
— Oq ■ ■ ■ ■ ^2r(fc-l) ~ [O.OO) 

Note that, with the Z^-symmetry constraints, 2rk Stokes matrices are reduced to fun- 
damental 2r Stokes matrices, {5'„}^^^, and also that the monodromy free condition is 
written as 

{st'^'^T-'Y = h. (3.54) 

4 The multi-cut boundary condition and solutions 

In the previous section, we developed the general framework of Stokes phenomena in the 
ODE systems which appear in the multi-cut matrix models. Mathematically, general 
solutions \E'(t; Q for these isomonodromy systems (or equivalently for the corresponding 
Douglas (string) equations) are parametrized by the Stokes multipliers with three basic 
constraints discussed in Section 1331 As is mentioned in Introduction, however, not all the 
solutions to these constraints can realize the critical points in the multi-cut matrix models. 
This consideration requires additional physical constraints on the Stokes multipliers. In 
this section, the first physical constraint is proposed, which we refer to as multi-cut 
boundary conditions. The second physical condition is proposed in Section [5l 
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4.1 Two different viewpoints about spectral curves 

Before we discuss the detail of the multi-cut boundary conditions, we first recall the set 
up of the multi-cut two-matrix models and the relationship between the Baker- Akhiezer 
function system (i.e. the ODE system) and the resolvent operator which defines the 
spectral curves. 

The definition of the multi-cut two-matrix models is given by the following matrix 
integral: 

dXdre~^*'^[^^(^)+^^(^)-^^l, (4.1) 



with the matrix contour C^"* of the following N x N fc-cut normal matrix, 

fc-i 

CjJ) = |x = f/diag(xi,X2,--- ,X7v)t/^; e f/(Ar), x^- G |J e^^^^t m}. (4.2) 

n=0 

The system of two-matrix models has the corresponding orthonormal polynomial system 
a„(x) = (^a;" + ■ ■ ■ ) , p^(y) = J=(y- + . . (4.3) 



with 



Sn,m = / rfxrfy e-^[^^(^)+^^(^^)-^^l anix) /3„(t/). (4.4) 



Here the contour C'-'^-* is given as 



k-l 



n=0 



= e IJe^-^tRl (4.5) 



an example of which is shown in Fig. [61 

The Baker- Akhiezer systems (or the ODE systems) appear as the double scaling limit 
of the orthonormal polynomials an{x) (or their dual polynomials /3n{y)), which is given 
as follows: 



a„ X = a 



"^/'^orth(C;t), (4.6) 



with the following scaling relations of a — )■ 0: 



X = ^0, ^ = exp(-to^) ^ 1, 

A^"' = ^str ^0, dr. = -a^'^g^trdt = -a^'^d ^ 0. (4.7) 

The continuous function \E'orth(^; C) is the scaling function of the orthonormal polynomials 
and satisfies the differential equations (14.81) and (14. 9p : 

C^orth(t;C) = i'(t;5)^orth(t;C), (4.8) 
^str^^orth(t; C) = Q{t] d) ^orth(t; C). (4.9) 
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(a) 



(b) 



Figure 6: Examples of contours C^*''-'. (a) is 6-cut contour C^-^^ and (b) is the 5-cut contour C^-^^ which 
is equal to the 10-cut contour C'^^^\ For reference, the position of cuts (zig-zag lines) around C — > oo is 
also denoted. 



This means that the orthonormal polynomial system is one of the solutions to the dif- 
ferential equations fl4.8p and (14 .9 1) , and eventually the ODE systems (I2.13P and (I2.14p . 
Consequently, the scaling function "^orthit; () is given by the canonical solutions \E'„(t; C) 



with some proper vector X^") = *(a;^' 



in) 



as 



^orth(t;C) = ^n(t;C)^^"\ 



[n 



0,1,--- 



(4.10) 



By taking into account the Stokes phenomena (I3.14p . these vectors X^"^ of various Stokes 
sectors Dr, are related as follows: 



(4.11) 



Note that the scaled orthonormal polynomials \l/orth(^;C) cire entire functions in G C 
because the original orthonormal polynomials are also entire functions. 

On the other hand, another important approach to solving the multi-cut matrix mod- 
els is the semi-classical approach with the resolvent operator Tl{x) of the matrix models. 



n{x) 



1 1 

— tr 

N X- 



X 



dz 

Qik) X 



(4.12) 



where p{z) is the density function of eigenvalues of the matrix X. An important fact 
about the resolvent is that this operator is a single valued function in x G C with a finite 
N and the cuts appearing in the large N limit are along the matrix-model contour C^^^ on 
the X space (as shown in Fig. [6]). These special cuts are called physical cuts. Interestingly, 
this resolvent operator is also related to the orthonormal polynomial solution \E'orth(^; C) 
in the following way 



(4.13) 
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^orth(t; C) ~ (det (x - X) ) ~ exp yx j dx'n{x') 
For the precise relations, see Appendix A in [35], for example. 
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with the scahng relation, x = u ^I'^dPI'^C^^ of Eq. (14 .Tp . This relation also indicate^ 

d 

n{x)^ lim (7strT^ln^orth(t;C)- (4.14) 

9str-s>o at, 

Therefore, the two different observables, the semi-classical resolvent 7l{x) Eq. fl4.12p and 
the semi-classical orthonormal polynomials ^orth(^; C) Eq. (I4.14p . provide two different 
viewpoints of spectral curves. 



4.2 The multi-cut boundary conditions 

As we carefully see the above two viewpoints of the spectral curve, one can notice that 
the realization of the position of physical cuts is not straightforward from the Baker- 
Akhierzer (or ODE) approach. The discontinuities of the scaling orthonormal polynomial 
of Eq. f l4.10p around C — )■ oo are Stokes lines and generally not distributed in the proper 
way expected in the semi-classical resolvent operator Eq. fl4.12p . This eventually means 
that not all the solutions to the ODE system Eq. (14. 8 p and Eq. (14. 9 p (therefore equiv- 
alently string equations) correspond to critical points of the multi-cut matrix models. 
The difference between these two viewpoints provides additional physical constraints not 
only on the vectors X^"^ but also on the Stokes multipliers si^ij which are identified as 
integration constants of the string equations. 

Next we formulate this physical constraint in the following way. Note that we here 
only care the leading behavior of ( ^ oo for the Stokes lines. 

Definition 10 (Multi-cut boundary condition) The Baker- Akhiezer (or ODE) sys- 
tems Eq. ( [^.g[ ) and Eq. ( [/^.g[ ) are said to satisfy the multi-cut boundary condition, if there 
exists a special solution \E'orth(^; C) which satisfies the following condition: 

• Stokes lines of the solution \l'orth(^; C) around C — ?• oo only exist along some special 
k angles ( ^ oo x e*^" : 

27m , , , , , 

X„ = Xo + -^ (n = 0,l,2,--- ,fc-l), (4.15) 

vjith a proper corresponding to each critical point. 

• Therefore, there exist an ordered set ofk indices, (ai, 02, ■ " " ? o-fc); and a set ofk non- 
zero vectors, (fi,f2,-"" yf^k), such that the asymptotic expansions of the solution 
^orth(^; C) '^n the angular domain D{xn,Xn+i) are given as 

^orth(t;C) ^ i;ne'^'""'(*^^) + --- , C^OOGL>(Xn,Xn+l), (4.16) 

asym 

and the expansions along the Stokes lines are given as the superposition: 

^oMt; C) ^ 6^'""'^*^^) + vn+, e^'^-^^H^) . . . ^ C ^ oo X e'^-^\ (4.17) 

asym 



Although ^orth(i;C) is a vector valued function, the behaviors of exponents are the same among 
the vector components. Therefore, it is understood by taking one particular element of the function 

^'orth(t;C)- 
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Here appears a special angle Xo which is determined by the critical points of the matrix 
models and is given as follows: 

— : Zfc-symmetric cases, and w^/^-rotated FSST cases / . , o\ 

k . (4.18) 

: Real-potential FSST cases 

This angle Xo comes from the scaling relation Eq. fl4.7l) . for detail discussion of which we 
should refer to [46j. An example of the boundary condition in the ( plane is shown in 

Fig.m 



/\\ 


(i;C) 




VWAW fi 


\ 

(i;C) ^ 





Figure 7: The multi-cut boundary condition in the 3-cut (1, 1) critical point. Although the general 
solutions to the Baker- Akhiezer function system can generally have "twelve cuts" , there are only three 
cuts in the scaling orthonormal-polynomial solution 5'orth(i; C)- 



Some comments are in order: 

• As in the definition of Stokes lines Eq. fl2.35p . we here only used the leading con- 
tributions of the exponents in — )■ oo. If we also take t — )■ oo, or equivalently 
if we just take Qstr — t- 0, on the other hand, we naturally encounter the general 
Stokes lines of Definition [51 Therefore, we interpret the general Stokes lines of the 
scaling orthonormal-polynomial solution as the non-perturbative definition of phys- 
ical cuts. In particular, this definition guarantees real eigenvalue- density function 
p{\)dX along the physical cuts: 

Re[7np{X)dX] = Re[d^^^\t; () - d^^^\t- ()] = 0, (4.19) 

where ( = ({\) is a local map from A G M to the generalized Stokes line GSLj C C. 
Note that this definition naturally justifies the curved physical cuts observed in 
[H] which appear when the matrix-model potentials are perturbed with complex 
coefficientso This consideration is further extended to off-shell backgrounds (or 
spectral curves ) in terms of the Riemann-Hilbert approach in Section El 

• In the p > 1 cases, the exponents ^p^^\t] have non-trivial cuts in the C plane, say 
f^^\t; C) ~ (^(9+1)/^. This p-th. root cut should be smeared by a proper supplement 

^^This consideration suggests that the position of physical cuts are not freely assigned and closely 
related to non-perturbative consistency with Stokes phenomena and therefore with D-instanton chemical 
potentials. 
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of exponents [S2]- This is also reviewed in Appendix |X1 Since we concentrate on 
the p = 1 cases in this paper, this point in the general fc-cut cases remains to be 
studied for future investigations. 

• As it will be clear in Section |4.2.2[ the set of indices (ai, 02, ■ • • , a^) in the multi-cut 
boundary condition (Definition [TOl) is generally given as 

an = j2r(n-l),fc = u + (u - l){r - 'j)mi, (4.20) 

with Theorem [3l In particular, the symmetric cases (7 = r) is given as a„ = n. 

Next we apply this boundary condition to concrete systems. Before devoting ourselves 
into general cases, however, we first consider the multi-cut boundary condition in the 
two-cut case, as a warm-up exercise for the general systems. 



4.2.1 The two-cut boundary condition 

Here we show how to solve the multi-cut boundary conditions in the two-cut (1,2) case. 
The orthonormal polynomial \l/orth(^; C) ^ Stokes sector Dn is generally given as a 
superposition of independent solutions, '^n\t;(): 



^orth(t;C) = ^n(t;C)X^"^ = xS")vi/W(t;C) + xr^i'^(t;C), C -> 00 g (4.21) 

However this assumption results in the 6-cut geometry of resolvent as shown in Fig. [8]-a, 
even though this system is called "two-cut". Therefore, one has to choose proper Stokes 
multipliers in order to satisfy the multi-cut boundary condition and therefore to obtain 
the geometry which only includes two cuts as shown in Fig. [8]-b. 
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Figure 8: The positions of cuts in the two-cut (1,2) ODE system, a) A general configuration of cuts 
for the general Stokes multipliers. There are 6 cuts, b) A configuration of cuts for the (1,2) critical 
point in the two-cut matrix models. The boxes indicate the regions Re(C) > and Re(C) < 0, in which 



the asymptotic expansion is given by ~ e 



(* 



position of cuts and the weak coupling infinity C 



1, 2). c) The profile of dominance depicted with the 
±00 e R. 



The multi-cut boundary condition is then given as follows: Since we wish to erase the 
cuts of orthonormal polynomial f l4.10p along the Stokes lines of 

^ = ±3, iy, (4.22) 
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we impose the following boundary condition on the vectors X^"); 

X">=(f), X>'>^('f]. X<^>=\%A. (4.23) 




where all the x^-"^ appearing here are non-zero. This can be also expressed in the dom- 
inance profile as in Fig. [8]-c. That is, if the Stokes sector Dn includes the following 
profile, 

[ mi I • • ■ I m/_i I mi | jhj+i | ■ • • | fhk-i \ fhk ] e D„, (4.24) 
then the boundary condition can be read as 

^oMt- C) = C), x^^l ^ 0. (4.25) 

i=i 

Since these vectors are related with the Stokes matrix (12.541) (with the Z2 symmetry 
condition (I2.59P ) as 

= 5„X("+^), X("+6) = X("), (4.26) 
one obtains the following conditions on the vectors X*^"^ and the Stokes multipliers: 



„(3)\ /^(4)\ /^(4)\ / ^(5) \ /^(5)\ / ^(0) 




which results in 



(4.27) 



/9' = 1, 7' = 1, l + /?7 = 0, 

) = xf ) = xl^) = xf ) = 7xf ^ 0, xf = = x« = xf = 7xf ) ^ 0. (4.28) 

Therefore, the solutions which are consistent with the Hermiticity condition ( I2.63P and 
with the monodromy free condition (12.661) are given as 

aeiR, /3 = -7 = ±1. (4.29) 

Consequently, the solution to the multi-cut boundary condition in the two-cut case has 
a real continuum parameter. However, as we will discuss in Section [5] with the Riemann- 
Hilbert approach, this parameter a causes "exponentially growing non-perturbative cor- 
rections" to the perturbative backgrounds (e.g. one-cut/two-cut spectral curves), except 
when a = 0. Therefore, the multi-cut boundary condition (with "the small instanton 
condition") completely fix the D-instanton chemical potentials as we advertised at the 
end of Section O 
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4.2.2 The multi-cut boundary-condition recursions (r = 2) 

From here, we solve the multi-cut boundary condition for an arbitrary number of cuts, 
k. In order to solve the constraints, we use the symmetric Stokes sectors (See Definition 



^orth(t;C) ^ ^2rl{t-X)X^^'''\ (^ooeD2rl, (/ = 0, 1, 2, 

asym 



(4.30) 



(sym) 



r ' Sq'^^'^ TK For sake of simplicity, however, we here 



Some of the results can be easily 



and its Stokes matrices, 5*2^; 
focus on the r = 2 cases, and therefore A; = 5, 7, 9, 
generalized to the general r cases. 

We first read the boundary condition in terms of the dominance profile: 

Proposition 1 (The multi-cut boundary condition on X*^"^) The multi- cut bound- 
ary condition in the general k-cut cases with r = 2 is given as 



The general k-cut cases 
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(4.31) 



Equivalently, the components of X^^"^ (r = 2 and k > 5) 



IS given as 



X. 



(An) 
n+i 



7^0 



X 



{An) 



k + l 



1,2, 



1,2,--- 



fc + 3 

4 . 

k + l 



(4.32) 



for n = 0, 1, 2, 



Eq. (gJip 



k — 1. The constraints on the Stokes matrices are then imposed by 



It is then convenient to introduce a new vector = (i/n,i)_i = r"X(4"), since the 
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Here fc = 3 is special because k < 2r = 4. This case is calculated separately in Appendix [ 
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above boundary condition becomes simpler: 



X 



(4n) 



7^0 



X 



„(4n) 



(4n) 



fc+3 



v 



y(4n) 



V 



(4.33) 



/ 



/ 



Note that the periodicity of index n follows: 

j^(4n) ^ j^{4(n+fc)) y(4n) ^ y(4(n+fc)) , . = y 



(4.34) 



In terms of the vector the constraints on the Stokes multipliers Eq. f l4.1ip are 

expressed as 



ji^{4n) _ ^(sym) j^(4(n+l)) 



^ = (5(«y-)r-i)FW«+i)). (4.35) 

Therefore, in terms of components, we obtain the following recursive relations for yn,i'- 

k 



I \ ^ (sym) 

yn,i — yn+l,i-l + 2_^^0,i,j yn+l,j-l, Vn+kJ — VnJ- 



(4.36) 



This is the central equations for the multi-cut boundary condition. After some tedious 
calculations, the multi-cut boundary condition turns out to be the following simple form: 



Theorem 5 (The multi-cut BC recursions) The recursion relation Eq. ^4-36 ) with 
the multi-cut boundary condition Eq. ^4-33 ) in the {k, r; 7) = (2m + 1, 2; 2) case is equiv- 
alent to the following two recursion equations for {yn,i}nez- 

LfJ L^J 

J^k[yn,l\ = Vn+m,! + ■Sl,m+2-j,l+j X 2/n+2i-l,l + S3^rn+3-j,l+j X ?/n+2j-2,l = 0, 

LfJ L^J 

Gk[yn,l] = -VnA + ^ So,fe+l-j,l+j X |/„+2j,l + ^ S2,k+2-j,l+j X |/„+2j-l,l = 0, (4.37) 

and linear expressions of the components {yn,i}i<f<k '^'^ terms of {yn,i}n&-' 

yn,i = yn,i{{yi,i}i&)- (4.38) 

Note that the coefficients in Eqs. i4-3'T^ are understood as modulo k, say S2,ij = S2,i+k,j- 
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The explicit expression for Eq. fl4.38p is a bit long and therefore shown in Appendix O 
with some examples. These recursive equations are the physical constraints which should 
be solved in addition to the basic constraints discussed in Section 13.31 In the general 
cases, the vectors Y^""^ in terms of {j/n.ijnez are denoted as 



(4.39) 



An important point is that all the Stokes multipliers si^ij in this expression are fine Stokes 
multipliers. Some detail derivation of this theorem can be found in jl02j . 

Finally we also make a comment on the boundary condition for general r (= 2, 3, ■ ■ ■ ). 
In terms of the dominance profile, they are expressed as 
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Here we define 



A^"'' = ji+2r{n-i),k = n+ [i/2\mi 

in) _ 



(4.40) 



(4.41) 



with Theorem [31 Note that A^ll^ = and aJ"^ 

sector D2r{n-i) includes the indices A'f^ of i = 3, 5, ■ 
profiles, then we impose 



A\ = n. Therefore if the Stokes 
, 2r - 1, 2r + 3, 2r + 5, ■ ■ ■ , in the 



X 



(2r(n-l)) 







3,5, 



and 



X. 



(2r(ji-l)) 
n+i 



7^0 



0,1, 



(4.42) 



The ending points of these series (about i) depend on how many segments -D2r(n-i) 
includes. This general classification could be tedious and we shall leave it for future 
study. 
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4.2.3 The complementary boundary conditions 

It is suggestive to sliow wliicli Stokes multipliers appear in the recursive equations 
Eqs. (14.37p . Here we show them by bold type in the profile of JT^^^™'' (i-^- Theorem 
H: 
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(4.43) 



Note that exactly a half of multipliers si^i^j <-> G ^7^^™'' picked up by the 

recursion. One may have the following question: Are there similar equations which pick 
up exactly another half of the multipliers ? This can be positively answered. In general, we 
can expect that there are (r — 1) similar equations, each of which picks up a different set 
of Stokes multipliers^ These equations come from complementary boundary conditions 
which are given by the multi-cut boundary condition (Definition [TUIl with different initial 
angles Xo: 



Xo 



TT 2'Ka 



1,2, 



!)• 



(4.44) 



The case of a = is the original multi-cut boundary condition Eq. fl4.18p . In fact, in 
the r = 2 case, the vectors in the recursion equation Eq. (I4.35P are replaced by the 
following F^^") 

/ VnA = \ 



y(4n) 



V 



(4.45) 



satisfying the complementary boundary condition of Eq. ( I4.44p (a = 1), and consequently 
we found the following different recursion equations with a = 1: 



Theorem 6 (The complementary BC recursion equations) The recursion relation 
Eq. lji4-36\ ) with the complementary boundary condition (a = 1 of Eq. ^.441 1) ^he 
{k,r;'y) = (2m + 1,2; 2) case is equivalent to the following two recursion equations for 

^^This anticipation is further explicitly shown in the fractional-superstring critical points {-f — r ~ 2) 
with arbitrary Poincare index r [103) . 
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{yn,m+2}ni=Z' 

LfJ 

J^k[yn,m+2] = ym+n,m+2 + ■53_/fc+2-j,m+2+j X 2/n+2j-l,m+2+ 

+ ■Sl,fe+2-j,m+l+j X yn+2j-2,m+2 = 0, (4.46) 

i=i 

LfJ 

Gk[yn,m+2] = ~yn,m+2 + ^ ^ S2^m+2~j,m+2+j X 2/n+2j,m+2+ 

i=i 

+ ■So,m+2-j,m+l+j X l/n+2j-l,m+2 = 0, (4.47) 

anc? linear expressions of the components {yn,i}i<f<,k terms 0/ {y„,m+2}nez-' 

yn,i = yn,i {{yi,m+2}i&z) ■ (4.48) 

Note that the coefficients in Eqs. (4^4^ ^^'^c understood as modulo k, say S2,i,j = S2^i+k,j- 



The explicit expression for Eq. fl4.48p is a bit long and therefore shown in Appendix O 
The vectors F*^"^ in terms of {yn,m+2}nei. also denoted by 

Y^''\{yi,m+2}iez) = {yn4{yi,m+2}iez))\ (4.49) 

In these recursion equations, exactly the complementary set of the multipliers si^ij -H- 
{j\i)i G Jk^2^^ are picked up. 

4.2.4 Useful reparametrization of the Stokes multipliers 

It is again suggestive to express these recursion equations as follows@ 

m 

[yn-m,l\ = yn,l + l/n-i,l = 0, 

i=l 



^^We point out the following interesting facts about these expressions. The algebraic equations defined 
by the recursion equations Eq. (|4.52p . 



m 

Gkiy) EE y-" g, [{y,- 1 y^},^^] = - (l + J] 0* y") = (4.50) 



n=l 



satisfy the following hermiticity relation: 

[My)]* = -y-'^Gkiy), if / = i. (4.5i) 

The same thing also happen for Eq. (|4.53|) . This therefore suggests that the solutions to the recursions 
{2/n,i}nez are given by fc-th roots of unity, y'' — 1. 

36 



lib 

Qk [yn,i\ = - (yn,i + 6** = 0; 

i=l 
m 

\_yn-m,m+2\ = 2/n,m+2 + di yn-i,m+2 = 

1=1 

m 

Gk \_yn,m+2\ = ~ ( Z/?i,m+2 + ^j* Z/n+i,m+2 



(4.52) 



0. 



(4.53) 



i=l 



The complex conjugation ^* (and ^*) comes from the hermiticity condition of Stokes 
multiphers Eq. ( I3.47p . It is also interesting to see the index n of the parameters 9n (and 
6n) in terms of the dominance profile: 
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(4.54) 



Here bold type is again the coefficients of the multi-cut BC recursions Eqs. f l4.37p . An 
important thing here is that these complementary boundary conditions are used to obtain 
explicit solutions of the Stokes multipliers (although they are not related to the physical 
boundary conditions). 

Finally, in order to write the explicit relation between the Stokes multipliers and the 
parameters 6'„, we introduce integers Li^ij G Z) which are defined as 



< Li^ij < k, Li^ij = (-1)' \i - j) mod k. 
In particular, we pick up the following set of indices (/; 

I - 1 



k = 4ko + l: Liu + 



G2Z+1; k = 4:ko + 3: Lii.+ 



(4.55) 



G 2Z, (4.56) 



and the relation is given as follows: 



Proposition 2 (The 6n parametrization) The fine Stokes multipliers si^ij are parametrized 
by k — 1 complex parameters {9n, 9n}nLi 

G J'ij"'^ satisfying Eq. ( [^.56] ) : Si^ij = 

and 

(iK) ^ ^i,!""^ satisfying Eq. ^.56\ ) : s^j 



(/ = 1,3) 



(4.57) 



^Lu,, (/ = 1,3) 
(/ = 0,2) 



(4.58) 



Therefore, this is a one to one correspondence up to the hermiticity condition Eq. ^.62^ . 
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4.3 Solutions in the general k-cut cases 

Before solving the boundary conditions, here we summarize the equations to be solved: 
After imposing the symmetry condition (13.39p . 

Zk symmetry: ^J^^ = T"' 5^ ^'"^ (/ = 1, 2, ■ ■ ■ , A; - 1) (4.59) 

the system becomes 

Multi-cut BC recursion: Y^^"'^ = (^f^^^^r-^) ^(^(n+i)) (4 gg) 

Monodromy free condition: (^Sq^"^^ T~^)'' = Ik, (4.61) 

Hermiticity condition: S'* = AT T^^A, (4.62) 

for n = 0,1,2, ■ ■ ■ , k — 1. In general, the patterns of solutions become complicated if we 
increase the number of cuts. However, we here show two kinds of special solutions which 
can be generalized to the cases with an arbitrary number of cuts. 

Before showing explicit solutions, we mention a key point of solving the above equa- 
tions. The main difficulty is from the monodromy free condition (I4.6ip . We here note 
that following fact@ 

Lemma 1 // the matrix Sq^^^^T~^ is diagonalizahle and its eigenvalues Xj are k-th roots 
of unity = 1, then the monodromy free condition Iji4-61\ ) is satisfied. The opposite is 
also true. 

Our strategy of finding solutions is now to show that Sq'^^^T^^ is diagonlizable. Below 
we list two types of explicit solutions: 

4.3.1 Discrete solutions and configurations of avalanches 

Theorem 7 (Discrete Solution) The following Stokes multipliers si^ij (written with 
0n,(^n of Proposition^^ are solutions to the multi-cut boundary condition in the sym- 
metric {p, q) = (1, 1) k-cut critical points (k = 2m + 1, 7 = r = 2^.- 

en = an{{-uj^^}f=,), en = an{{-oj^^}f=,), {n = l,2,--- ,m) (4.63) 
with the symmetric polynomials an among {xi}fLi of degree n@ 

Cni{xi}i=i) = ^ ^ ^^jjXjj ■ ■ ■ ajj^ , (4.64) 

l<il<i2<-<in<N 

if and only if the integers (ni, ^2, ■ ■ ■ , n^; ni, n2, ■ ■ ■ , n^) of Eq. 1^4 -S^S^ satisfy 
rii ^ rij. Hi ^ rij mod k {i ^ j). 



— "^^{nj + Uj) ^ Hi, Hi mod k (i = 1, 2, ■ ■ ■ , m). 



(4.65) 



^^The opposite is non-trivial but can be shown by using Jordan normal form. 

^^Here cr„ stands for the symmetric polynomials. Do not be confused with the Pauli matrices. 
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Comments on the conditions Eqs. ( I4.65P are following: 

• In this solution, one can find the (k—l) explicit eigenvectors of the matrix S'g'^'^^^r"^: 



(4.66) 



with the vectors of the BC recursion equations, Eq. f l4.39p and Eq. fl4.49p : 

y,, = y W ^ u;"""^}„ez] , yj = [{yn,m+2 ^ cu^'^^nez] . (4.67) 

They are distinct only when ^ rij {i ^ j) and ^ rij {i 7^ j) with modulo k. 



Noting that det S^ 



{sym)p_i 





1, one concludes that the eigenvalue of the remaining 



eigenvector is given by uj "° with hq = — + ^7 ) • This eigenvector becomes 

distinct only when Ui ^ uq ^ nj = 1, 2, ■ ■ ■ , m)o 

With these considerations, one can prove the above theorem. Below we show a graphical 
expression of the conditions Eq. fl4.65p in terms of Young diagram. 

1. The following transformation is an automorphism among the solutions to the con- 
ditions (14:651) : 



{ni,n2, 



, n. 







(ni + 1, n2 + 1, ■ ■ ■ , + 1; ^1 + 1, n2 + 1, 



+ 1), (4.69) 



which also maps Hq (= — X]j=i(% + ^i)) Uq nQ + 1. 



2. By choosing the following representative of the solutions as no = 0, and by properly 
choosing the ordering of the indices, one can rewrite the conditions f l4.65p as 



J](n,+n,) =0, 



\ < rii < n2 < ■ ■ ■ < ri|^fej < — 1, 
1 <ni <n2 < ■ ■ ■ < n^fej < k — 1 



(4.70) 



3. Therefore, these indices can be expressed in terms of Young diagram. Here is an 
example (a solution in the 11-cut case). 



(4.71) 




(ni, n2, ^3, ^4, ns) = (1, 2, 4, 6, 9) <^ 
That is, the i-th row from the bottom has Ui sky boxes in the diagram. 



^® As a side remark, here we show the eigenvector J^o of the eigenvalue 77 ^ = lj "° when rij = rij [j = 
1,2,--- ,m: 

= ,ryL¥J,o,... ,0,(-l)V/',--- , (-1) V^'+L^J , 0, • • • , o) *. (4.68) 



IV 



Note that ah the components of Vo in the region II and region IV vanish (See Eqs. (jC.l|) for definition 
of the regions). 
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We also draw mxk total boxes for later convenience. In particular, the upper-left Young 
diagram (written with ^) is referred to as sky and the lower-right Young diagram (written 
with □) is as snow. The pair {nj;nj)j is denoted as 

(ni, n2, ns, n4, n^; ni, n2, na, n^, n^) = (1,2, 4, 6, 9; 3, 5, 7, 8, 10) 




(4.72) 



Therefore, the graphical meaning of the conditions Eqs. fl4.65p is following: 

• The number of the boxes EH (amount of snow) is always a multiple of fc, and the 
following configurations are allowed solutions in the 7-cut case: 




(4.73) 



• Neither rii and (z = 1, 2, ■ ■ ■ , m) can be or A;, therefore the following configu- 
rations are not allowed: 




(4.74) 



The solutions cannot have vertical cliffs, therefore the following configurations are 
not allowed: 



forbidden: 







(4.75) 



One of the ways to exhaust the solutions is first to take the most steepest configurations, 
and then to consider possible ways for snow to slide on the surface with satisfying the 
condition ( 14. 70 p . for example: 
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(4.76) 



We refer to these configulations of Young diagram as avalanches. Therefore, we conclude: 
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Proposition 3 (Avalanches) The discrete solutions to the non-perturbative completion 
are labeled by configurations of avalanches in terms of Young diagram. 

Note that one can also move some snow on the one side to the other side. It is also worth 
mentioning the following two transformations: 

• The first transformation is called dual, 

(ni,ra2, ■ ■ ■ ,nm;ni,n2, ■ ■ ■ ,nra) 

[k - ni,k - n,,- ■ ■ ,k - n^;k -ni,k -n,,- ■ ■ ,k -n^), (4.77) 

In the terminology of Young diagram, the dual transformation f l4.77p exchanges the 
sky and snow of the left and right Young diagrams simultaneously. In particular, 
the following diagram shows the action of the dual transformation on the left Young 
diagram: 




(4.78) 



• The following is called reflection, which exchanges the snows of the left and right 
Young diagrams: 

{ni,n2,- ■ ■ ,nm]ni,n2,- ■ ■ ,nm) (^i, ^2, ■ ■ ■ , n^; ni, ^2, ■ ■ ■ , n^) . (4.79) 

Note that these two transformations also automorphisms which fix the condition (14.701) . 

4.3.2 Continuum solutions 

Theorem 8 (Continuum Solution) The following Stokes multipliers si^ij (written with 
On, On of Proposition\^ are solutions to the multi-cut boundary condition in the sym- 
metric {p, q) = (1, 1) k-cut critical points (k = 2m + 1, 7 = r = 2j.- 

On = crni{-u^^}]L,), On = Sn{{0,}]l,)+d*^_^^,0*^, {u = 1, 2, ■ ■ ■ , m) , (4.80) 
with the polynomial Sn{x) which are defined by the following recursion relation: 

n 

^n{{xj}jez) = - ^Xi5„_i({a;j}jgz), So{{xj}jez) = 1, (4.81) 

i=l 

if and only if the integers (ni, n2, ■ ■ ■ ,n^k^) satisfy rii ^ rij mod k {i ^ j). 

A derivation of this solution is shown in Appendix |Dl The concrete expression of 
Eq. (I4.80p (and therefore the polynomials Sn{x)) is given as 



Oi = 


{-Oi) + 0*n^O*n^, 


02 = 


(-6*2 + Ol) + Oln^^O*^, 


03 = 


(-6*3 + 2O1O2 - of) + O*^_20*m^ 


0,= 


[-0, + Ol - W\02 + 20,0^ + of) + Ol^_^0*^, 


0,= 


(-^5 + OS, - 'iOlO-^ + O2O3 + ^0\02 - 30,01 - ^1) + 
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It is worth mentioning the relation to the Schur polynomials -P„, ({xjjjg^): 

Sn{{xj}jez) = Pn{{yj}j&z), Xn = Pn{{-yj}jez), (4.83) 
where the Schur polynomials Pn{{xj}jez) are defined as 

oo oo 

J2^''Pn{{xj},e7) = exp[^z'^a;„]. (4.84) 

n=0 n=l 

Note that these solutions includes m (= [|J) real parameters. Sometimes, eigenvalues of 

the matrix Sq^^^ of the discrete solutions are distinct. In this case, such a discrete 
solution is a special case of the continuum solution. However generally these solutions 
do not include the discrete solutions in Section 14.3. ![ since the discrete solutions gen- 
erally include degeneracy of eigenvalues which cannot be resolved by these continuum 
parameters. 



5 Stability of perturbative backgrounds 

In this section, we briefly review the Riemann-Hilbert approach and the Deift-Zhou 
method [8H - I86] . and also discuss its physical interpretations in non-critical string the- 
ory. In particular, we argue that this procedure implies an additional physical requirement 
about stability of classical (or perturbative) backgrounds. We will see that this constraint 
results in the proper Stokes multipliers expected in the two-cut (1,2) critical point. Clas- 
sical background here means the spectral curves which appear as semi-classical (large N) 
solutions of matrix models. 

The role of the Riemann-Hilbert approach is to obtain the t dependence of physical 
amplitudes (for example, asymptotic expansion in t) by using an integration expression 
which can be derived from the ODE system in The review article [91] contains useful 
references of the Riemann-Hilbert approach. 

Roughly speaking, in the Riemann-Hilbert approach, we flrst discard the analytic 
continuity of the canonical solutions f l2.38p and keep the form of asymptotic expansion 
( I2.25P in the complex plane C. In practice, we introduce some Stokes sectors (here we 
consider flne Stokes sectors) Dn and canonical solutions on them, \E'n(t; C)- As it has been 
reviewed in Section 121 these canonical solutions have the same asymptotic expansion in 
each Stokes sector f l2.38p and the difference of these canonical solutions is expressed 
by Stokes matrices fl2.39p . Therefore, inside the intersection of two Stokes sectors, we 
introduce a semi-inflnite straight line from the origin, /C„, 

)Cn = {C = ue'^-,ueR+} with ^Xne[0,27r) s.t. /C„ C D„ n (5.1) 

and deflne the following new function \E'rh('^; C) which is analytic in ^ G C\ K, with 

^RH(t;C) = ^n(t;C) C e D{Xn-l,Xn), (n = i,2,---), (5.2) 
which has the following uniform asymptotic expansion in C G C \ [j^JCn- 

^RH(t;C) ^ $asym(t;C) = ?(t;C)e^(*^^\ C^oogC\M/C„. (5.3) 
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The lines /C is referred to as discontinuity lines, and examples are shown in Fig. [HI Note 
that the function \I/rh(^; C) has enough information to recover all the canonical solutions 
simply by analytically continuing the argument (. 




Figure 9: These are examples in the 3-cut (1, 1) critical point, a) The coarse Stokes sectors (shadowed 
domains) and the discontinuity lines /C (dashed lines). Basically, any lines in the intersections D^n H 
^sin+i) E^re allowed, b) The discontinuity lines K. (dashed lines) with respect to the fine Stokes sectors. 
They are related to the lines in (a) by continuous deformations which do not cross any divergence in the 
Riemann-Hilbert integral (|5.10p . 

An essence of the Deift-Zhou method for the Riemann-Hilbert problem [86] is in- 
troduction of the following k x k function g{t]Q which we shall call (off-shell) string 
background: 

g{t;0 = diag{g'^'\t;Or-- ,g^'\t;0), 

r CO _ 

With /)(t;C) = $^t«C + ^^lnC + E-^n^C-^ (5.4) 

n=l n=0 

If one focuses on the aspect of algebraic curves, the function g(t; () is referred to as (off- 
shell) background spectral curve. We then obtain the following setting of the Riemann- 
Hilbert problem: 

Lemma 2 (Setting of the Riemann-Hilbert problem) There exists the set of pa- 
rameters tn"* (i = 1, 2, ■ ■ ■ ,k; n = 1, 2, • • • , r) which satisfies 

Zit; ^ ^Ruit; C)e-^(*^^) ^4, (( ^ oo G C \ |J ^n) • (5.5) 

n 

The k X k matrix function Z{t; () then satisfies the following discontinuity relation: 

Z+{t- = Z_{t; Gn{t; C), G„(t; C) = e^(*^«S'„e-^(*^^), along ( G /C„, (5.6) 

where n = 1, 2, ■ ■ ■ and we define Z±{t; () = lima^o Z(t; ( ± ae) with a vector e which 
directs towards the left hand side of the line /C„. 
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In general, the parameters t 



1,2,--- ,k; n 



1, 2, ■ ■ ■ , r) are the integrable de- 



formations of the fc-component KP hierarchy [96j. These are then given by the Lax 
equations: 

d 



5'str- 



dt 



This information is understood as given information of the system and non-normahzable 
string moduh space which should not be minimized by the string dynamics [81]. Note 
that the Stokes matrices are invariants of these integrable deformation: 

dS„ 



dt 



(i) 



0, 



k; n = 1,2,--- ,r; m = 1,2,---), 



{51 



and therefore the multipliers are integration constants (initial conditions) of these defor- 
mations. In this sense, they are also understood as non-normalizable string moduli space 
of the dynamics in the strong-coupling region of string theory. 

Since the Stokes multipliers are integration constants of the system, we can uniquely 
obtain all the information by identifying the deformation parameters tn^ {i = 1,2, - ■ ■ ,k; n 
1, 2, ■ ■ ■ , r) and the Stokes multipliers. The fact is given in the form of the following the- 
orem: 

Theorem 9 (The Riemann-Hilbert problem (see |91| )) For a given analytic func- 
tion G{t; () on the discontinuity line ( & K, = [J^ /C„, 

G{t; C) = Gn{t; C) = e^(*^^)S„e-5(*^^) ( E )Cn (n = 1, 2, ■ ■ ■ ), (5.9) 

there exists a unique holomorphic function Z{t; Q which satisfies Eqs. / I5. 5\) and / (5. 6]) . 
and is given as 

d\ p(A)(G(A)-4) 



Z{t;0 = h+ [ 
Jk. 

= 4 + E 



27ri 
s 



A-C 

27ri A-C 



(5.10) 



wtthp{C) = Z_{0 on Cg/C = U„/C„. 

By using the Riemann-Hilbert solution f l5.10p . one can obtain the canonical solutions to 
the ODE system (defined in (15. 3p ) as a function of t: 



^R,H(t;C) = ^(t;C)e 



9{t;C) 



asym 



asym 



{t-C) = Y{t-C) e 



C^oogC\/C. (5.11) 



Note that the "density function p(C)" is given by Z{t; () itself, and then the function p{Q 
satisfies the following integral equation: 



p(C) 



dX p{X){G{X)-h) 



Ce/C. 



(5.12) 



27ii A — C + e 

Therefore, one can recursively solve it and the solution is given as the following infinite 
sum of integrals: 



n=l i=l L""^ J j=2 



Xj — Xj—i -\- 6 



G{Xi) - h 
Ai-C 



(5.13) 
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with the assumption that 

/ ^(G(A)-4), (5.14) 

is sufficiently small. Note that we use the following multiplication rule of matrices: 
YYj=i Aj = AnAn_i ■ ■ ■ Ai. In terms of componets, this is expressed as 

Z{t; = 4 + 2^ Sn,ijE, 



2m Ai - C 



+ 



-r;^ f dXi f (i\^(.9^^Kt-M)-g^'Ht-M)+9'^'Ht-M)-a'^'Ht-M) 

+ 2^ Sr,,,i,lS^,,ljE,, j — — ^ +.... 

ni,n2,i,j,l •^'^"1 "^'^"2 ^ ^ ^ 7V i sy 

(5.15) 

This expression is formally convergent if the subsequent integral are small enough. In this 
case, one can evaluate the leading contribution by truncating higher terms (the so-called 
Born approximation). It is worth mentioning that this integral is quite similar to the 
D-instanton operator formalism in the free-fermion formulation [3T|[32] by interpreting 
g'^'^^t; () as the free boson operator v^o''(C) in the system. 

An important point here is that, in the Riemann-Hilbert approach, the string back- 
ground g{t] () is arbitrary except for the parameters tn\i = 1,2, ■ ■ ■ ,k] n = 1,2, ■ ■ ■ , r), 
and then generally is different from the semi-classical resolvent amplitudes ip{t] () of 
Eq. f l2.32p which is obtained as a solution to the equation of motion (or loop equations) 
in the large N limit of the matrix models. As one can see in Theorem [9l the role of 
the string background g{t; () is a reference background in the Riemann-Hilbert problem. 
Therefore, from the string-theory viewpoints, the string backgrounds g{t] () are generally 
understood as off-shell backgrounds of string theory and in this sense the Riemann-Hilbert 
approach realizes an off-shell background independent formulation of string theory. 

In order to understand g(t;() as off-shell backgrounds of string theory, it is worth 
mentioning the interpretation of the position of cuts. Taking into account the consider- 
ation given around Eq. f l4.19p . we can define the cuts on the off-shell background as a 
combination of general Stokes lines: 

Re(^«(t;C)-^7(^ni;C)) =0, (5.16) 

which is obtained by an analytic deformation of the matrix contour u^^'^C^''^ (so that it 
realizes the multi-cut boundary condition around ( — t- oo). Note that this consideration 
is possible after imposing proper Stokes phenomena which solve the multi-cut boundary 
condition, as it is carried out in Section SJ 

This viewpoint also provides the following consideration: If one chooses g{t; () as a 
semi-classical resolvent function ip(t; (), then the evaluation of Eqs. (15. lip and fl5.10p in 

5'str 0, 



^Mf = Z{t; C) e^^'-'^^ = [4 + ■ ■ ■ J e^^'-''^\ (5.17) 

is a calculation of quantum corrections from the background spectral curve g{t; Q which 
is given by the semi-classical resolvent. Therefore, if the resolvent background is a stable 
vacuum of this system, the non-perturbative corrections should be exponentially small. 
This is the additional constraint for the Stokes multipliers and is referred to as small- 
instanton condition. 
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5.1 The small- instant on condition for the 2-cut critical point 



Here we consider the small-instanton condition in the 2-cut (1,2) critical point. Mathe- 
matically, the Riemann-Hilbert problem in this case has been evaluated in [85 | [86 |[92] - [95] 
in the larger classes of Stokes multipliers (See the review [9T]). In particular, according to 
the Deift-Zhou procedure [HE], one first deforms the discontinuity lines /C to anti-Stokes 
lines. The concept of anti-Stokes lines depends on saddle points of the string background 



saddle points C* 



d 



0, 



c=c* 



ihJ = 1,2,--- ,/i;). 

(5.18) 

Definition 11 (Anti-Stokes lines) Anti-Stokes lines ASL^"'' are defined for each pair 
of {id) as 



ASL 



(n) 



CGC;Im g^''^\t-0 



Im 



(5.19) 



where Q" is a saddle point of the function g^'^'^\t; () = g^'''\t; Q — 9^-'\t] Q. 

In the procedure of the Deift-Zhou method, one can choose the string background g{t\ C), 
however, we know that the 2-cut (1,2) critical point has two phases with respect to the 
sign of t cosmological constant ^2]. Therefore, we choose the string background according 
to the actual phase appearing in the two-cut matrix mo del 



9{t\ C) = 



3C» + «C + 



0-3 
0-3 



: two-cut phase (t > 0) 
: one-cut phase {t < 0) 



(5.20) 



Since we know that these curves are realized in the critical point as its stable vacua, these 
perturbative vacua should satisfy the small-instanton condition. Below we consider each 
case separately. We skip the calculation which is the same as that in \91\ . 



The two-cut phase (t > 0) There are three saddle points of the function g^^''^\t; Q = 
9^'\t;0-9^'Kt;Cy. 

C = C;5: CS = 0, Cif = ±zV2t, (5.21) 
and the values of the function at these saddle points are 

9^'^'\t;Ci'l) = l{2tY/\ g^'''\t; C[f^) = 0. (5.22) 

Note the saddle-point value of the function g^'^'^\t; () = —g^^''^\t; Q. They are understood 
as instanton actions for the saddle points. The deformation of discrete lines /C to the DZ 
curves is given in Fig. [TOl 

^^Note that we are here imposing a physical requirement, by taking into account the Deift-Zhou 
method [86]. In the Deift-Zhou procedure, one considers an arbitrary Stokes multiphers, and the function 
g{t; C) is a function which we choose so that there is no divergence in the RH calculation. In this way, we 
can obtain the asymptotic form in t for these arbitrary Stokes multipliers. In this section, on the other 
hand, we impose a physical constraint in which the physical background g{t; C) obtained from the matrix 
models is stable perturbative background with small non-perturbative effects. Therefore, this constraint 
picks up the special and physical Stokes multipliers. 
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(a) (b) (c) 

Figure 10: The discontinuity lines and the DZ curves in the two-cut (1, 2) critical point of the two-cut 
phase, a) The discontinuity lines K,. There are two kinds of lines: the one kind is the lines K,2n+i on 
which the integral (|5.15p only includes the contributions from the exponent e^' The other kind is 

the lines K,2n on which the integral (I5.15|) only includes the contributions from the exponent e^'^'^''^^'. b) 
The DZ curves which are obtained from analytic deformation of the original lines /C. A large D-instanton 
effect appears around the origin on the line /C3. Therefore, we require a = so that this large instanton 
vanishes, c) The resulting DZ lines with a = 0. Two lines along the real axes /Co± come from the 
Stokes matrices on the lines ICjj and ICd- Saddle point approximation on each line gives ZZ branes in the 
Liouville theory, however contributions from these lines are the same and canceled by the Z2 symmetry. 



On the DZ curves, we then evaluate the integral flS.lSp at saddle points [9T]. The 
small-instanton condition becomes relevant when the saddle point ^1*2 = of g^^''^\t;C) 
contributes in the Riemann-Hilbert integral flS.lSp . This happens in the integral on the 
curve /C3. The relevant part is given as 

Z{t;0 = aE^,, / ^^^ + .... (5.23) 

The parameter a is the Stokes multipliers of this system fl2.59p . Therefore, the small- 
instanton condition requires 

a = So = S3 = 0, (5.24) 

otherwise this perturbative vacuum (15.201) breaks down (or decays into some stable vac- 
uum) by the large non-perturbative effects. Consequently, the solutions to the non- 
perturbative completion are finally fixed to be 

a = 0, /3 = ±1 = -7, (5.25) 

which is known as the Hastings-McLeod solution in the Painleve II equation [89j. As it 
has been calculated in ^9J, the final result is given aS 

f{t) = -2f3V2i+--- with /3 = ±1, (5.26) 



'^'^In this calculation, we use the local Riemann-Hilbert problems. Since the evaluation of the Riemann- 
Hilbert problem is not our purpose, we here skip the calculation. See the review |91j . An intuitive reason 
for vanishing the D-instanton effects (or physical interpretation of the mathematical result) is cancellation 
due to the Z2 symmetry of the system. For example, if one introduces the formal monodromy (as 
mentioned around (|2.64p . i.e. adding DO-brane charges in the background) then the instanton effect from 
the origin C = appears. 
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especially, the instanton effect which comes from a single ZZ-brane at the origin C = 
vanishes in this phase. 



The one-cut phase (t < 0) There are two saddle points of the function g^^''^\t; Q = 

C = CS : Cif = ±v^, (5.27) 
and the values of the function at these saddle points are 

9^'^'\t;Cif^) = T^{-t)y\ (5.28) 

The deformation of discrete lines K, to the DZ curves is given in Fig. [11] Note that 
existence of this phase also requires the same constraint a = 0. By taking into account 
the solution to the non-perturbative completion f l5.25l) . the Riemann-Hilbert integral 
(I5.15P becomes the following simple contour integrals: 

Jk^ 2 2m A - C Jk2 1 27rz A - C 



= h^ "^W — 1= , i\ — ■;= ■;= 

therefore the asymptotic expression of f{t) is given as 



e-K-*)^^^ + ■ ■ ■ , (5.29) 



f{t) = ^ e--^i-tf^' + --- with /3 = ±1. (5.30) 

V27rv/^ 



See Eq. fl2.48p . It is worth mentioning that a similar expression was found in the 2-cut 
(1,2) critical points [71] which comes from an explicit expression of fermion state within 
the free- fermion formulation [25l|3Tl[32l|69] , although the expression there is given by an 
infinite sum of super-matrix integrals. 



5.2 The small-instanton condition for the /c-cut critical points 

Here we consider the small-instanton constraint in the /c-cut (1, 1) critical points. Since 
we here focus on the additional constraint, we only study the saddle point actions for 
the semi-classical string background and evaluation of the Riemann-Hilbert integrals is 
remained for future investigation. The classical backgrounds in these cases are calculated 
in [l6j and given in terms of parameter z as 

^7(t;C) = diag(^7«(t;C),---,^7^'Hi;C)), 9^'\t-X) = f \{z)dx{z), 

with x{z)=t{J{z-c)\z~hf~\ y{z)=t^J {z-cf'\z-h)\ (5.31) 

with = cl+h {k — l). The index Z (= 0, 1, 2, ■ ■ ■ , k — 1) labels generally different solutions. 
The classical background g{t; Q is then expressed as 

9^'\t- C) = (t; u-^^-^K) , 9^'^^ {t; x) = i (zix))' -ic + b) z{x). (5.32) 
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Figure 11: The discontinuity lines and the DZ curves in the two-cut (1, 2) critical point of the one-cut 
phase, a) The discontinuity lines /C which is the same as two-cut phase, b) The DZ curves which are 
obtained from analytic deformation of the original lines /C. A large D-instanton effect appears around the 
saddle point ( = +\/t on the line /Co, and around the saddle point = —Vi on the line /C3. Therefore, we 
require a = so that these large instantons vanishes, c) The resulting DZ lines with a — 0. By taking 
into account the sign of the Stokes multipliers, one observes that the integral (|5.15p along connected 
lines /C2 and /C4 (and also /Ci and /C5 in the same way) can be considered as an integral on the single 
contour. Saddle point approximation on each line gives ZZ branes in the Liouville theory of the one-cut 
phase. 



Here z{x) is the inverse of the function x{z) in Eq. fl5.3ip . The saddle points for 
g(''^\t; C) = g^'\t; () - () are given as 

^J''^\t- 0=0 ^ u^-'x{z) = u^-'xiz'), u-^^-M^) = u-^^-'k{^'), (5.33) 

and then this can be solved as 

2cos(xS'/2) )• = -I 2cos(x!"'/2) J' ' ' 

with xtj = 2vr^^^i3^^, (n = 1, 2, ■ ■ ■ ). Substituting these values in Eq. f l5.32p . we obtain 
the saddle point action: 



'^''-iz^"')')-(b + c)(z^-z^-^ 



tan(^) G iR. (5.35) 

This means that the saddle point action always contributes order 0{g^^^^) and then iden- 
tified as perturbative corrections (not as instantons). Therefore, there is no additional 
(small-instanton) constraints on the solutions obtained in Section |H 



6 Conclusion and discussions 

In this paper, we give concrete solutions to the non-perturbative completion in the fc-cut 
two-matrix models by a quantitative study of Stokes phenomena. The non-perturbative 
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completion problem consists of the multi-cut boundary condition for the orthonormal 
polynomial systems and the non-perturbative stability condition for the semi-classical 
spectral curves in the large N limit. By carrying out these procedures, we demonstrated 
two classes of solutions, which are referred to as discrete and continuum solutions. Inter- 
estingly, the solutions possess kind of "charges" in terms of Young diagram representation. 

We note that the continuum solutions to the non-perturbative completion still include 
continuous free parameters, although the two-cut cases have been completely fixed. It 
is conceivable that we might need to rely on further independent physical arguments to 
reduce these degrees of freedom, here we would like to interpret these free parameters 
as physical moduli parameters in the non-perturbative region of the string theory. Since 
the strong-coupling dual theory of the multi-cut matrix models seems to be non-critical 
M theory [44j, these continuous parameters would correspond to the non-perturbative 
(non-normalizable) moduli space of M theory, A^M°theory™ which is a distinct parameter 

space from the string-theory moduli space, A^itring '^"'^^ ^^"^ -^stxin™ Below we provide 
a list of issues which deserve further exploration. 

• In this paper, we have solved Stokes phenomena in Z^-symmetric critical points. It 
is also interesting to consider similar program in the fractional-superstring critical 
points [13] . In particular,we would like to see the emergence of the non-critical M 
theory from the k ^ oo limit 

• Our procedure is directly related to Riemann-Hilbert calculus. It is useful to ex- 
amine higher order instanton sectors and generalize the results in [78j. 

• In this paper, we focus on the cases with p = 1 and small q. In order to extend this 
procedure to the general q cases, one should resolve several complexities as shown 
in Eq. f l4.40p . It is of great interest to obtain the Stokes multipliers in higher {p, q) 
critical points. In particular, evaluation in the bosonic cases would clarify the issue 
raised in [61] . Also we have to take into account the smoothing of the cuts as shown 
in [52] (also see Appendix lAj) . 

• It is interesting to investigate whether the Riemann-Hilbert representation can be 
written in language of matrix models? This resembles the supermatrix models 
[7T] which appear by evaluating tau-function in terms of free fermions. Also it is 
interesting to compare it with Kontsevich type matrix models [105 J and also with 
the non-perturbative topological string-theory block recently proposed in [106] . 

• The Riemann-Hilbert representation is a background independent formulation, which 
allows us to introduce general off-shell background in string theory. Therefore, it is 
interesting to study physics in off-shell backgrounds and general concept of back- 
ground independence in matrix models/string theory. 

• In the multi-cut matrix models, there are two kinds of perturbative string vacua |44j : 
One is perturbatively isolated sectors (perturbative superselection sectors) which 
are decoupled with other sectors in all-order perturbation theory. This phenomenon 
is an origin of the extra-dimension in M theory. The other is perturbative vacua 
in the string-theory moduli space. For survey for the second vacua, the Riemann- 
Hilbert representation is even more powerful, since the off-shell moduli space is 
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understand as the space of off-shell string-theory backgrounds. Furthermore, the 
Zfc symmetric critical points in the multi-cut matrix models have several pertur- 
bative vacua which satisfy loop equations. Therefore, it is interesting to study 
non-perturbative string-theory landscape from the Riemann-Hilbert approach. In 
particular, it might be possible to identify which observables are suitable for a 
discription of a potential picture in the moduli space. 

• We obtained several solutions to Stokes phenomena which are characterized by 
several charges carried by Young diagrams. What is the physical meaning of our 
solutions? Any relation to W-symmetry or WZNW? 

• Our solutions are natural generalizations of the Hastings-McLeod solution in the 
Painleve II equation. The Hastings-McLeod solution is known to have several 
special features, for instance analyticity of the solution (See also [91j). Therefore, 
it is mathematically interesting to understand the analyticity of the solutions in t 
and to identify the standing point of our solutions in general solutions of the string 
equations. 

• As is well-known, the integrable deformations in the usual integrable system corre- 
spond to the moduli space of worldsheet conformal field theory. On the other hand, 
non-trivial deformations of our solutions can be interpreted as non-perturbative in- 
tegrable deformations in physical solutions of string equations. Therefore, these 
deformations are related to the moduli space of the dynamical degree of freedom 
in the strong coupling region, i.e. degree of freedom in non-critical M theory. It 
is interesting if there is a comprehensive understanding of these non-perturbative 
integrable deformations. 
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A Stokes phenomenon in the Airy function 

The non-perturbative relations between the resolvent and the orthonormal polynomials 
are first studied in [52] in the (2, 1) critical point of bosonic minimal string. Since this 
study also uncovers another aspect of cuts in the resolvent curves for the cases of p > 2, 
we here briefly review the results and summarize the key points. 
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In the bosonic (2, 1) critical point, the orthonormal polynomials satisfy the following 
differential equation: 



C^orth(t; C) = {d' + u{t)) vi>„rth(t; C), (A.i) 

(7str^^orth(t; C) = 5^orth(t; 0- (A.2) 

By taking into account the definition d = Qstrdt, one can show that the orthonormal 
polynomial is given as Airy function: 

= [gL ^°rth(t; C), *orth(t; C) = Ai(C + 1). (A.3) 

Here we have concluded u{t) = —t by imposing the integrability condition of flA.ip 
and ( ]A.2p . and also have chosen the damping solution (Airy function) as the physical 



solution [52] : 

^orth(t; C) ^ 0, C ^ OO. (A.4) 

As it is well-known, the asymptotic behavior of the orthonormal polynomial \l/orth(^; C) 
(i.e. the Airy function) around the real axes, C, — t- ±oo, is given as 



in ^ — )■ oo with the angle, | arg(C)| < tt, and 

^orth(t;C) - I. , e ^sstr^^ +«e39st.^^ +••-, (A.6) 



in — )■ e'^* X cx) with the angle, |arg(— (^)| < 27r/3. Note that both two expressions 
in the intersections, vr/S < |arg(^)| < vr, have common asymptotic expansions, and 
therefore, the appearance/disappearance of different exponents in different asymptotic 
regions is understood as the Stokes phenomenon. As a consequence, the resolvent in the 
weak coupling limit (^str — ?■ is smooth in ^ with arg((^) < tt, and the discontinuity only 
appears along ^ G (—00, — t), that is. 



lim 

e->-±0 



lim ^orth(i;C + «e) 

3str->0 



~ e^^M''^'"\ C e (-00, -t). (A.7) 



An important point in [52] is that the resolvent curve itself has a cut around ( ^ 00. 
However the explicit cuts are smeared by the superposition of the exponents e*^''"''*^^'^^ and 
g-(C+t)^'^^_ Note that the solution corresponding to matrix models can be chosen by the 
single condition Eq. ( ]A.4p . This is due to the simplicity of Airy system. In more general 



system, however, one needs the multi-cut boundary condition proposed in this paper. 

B Lax operators in the multi-cut matrix models 

Here we summarize the Lax operators used in this paper. 
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B.l The Zfc symmetric (1, 1) critical points 

This class of critical points are characterized by the following Lax operators: 

P{t-d) = Td + H{t), 

Q{t- d) = (r-2(t; d) Pit; d))^ - fi{T-\t; d))^ 
= T-^d~T-^HT-^ - lj,T-\ 



(B.l) 



Note that the symmetry requires 



H{t) 



/O * 

* 



* 



(B.2) 



and the Lax operator r(t; d) is defined as 

oo 

^(^;9) = ^ + ^5„(t)9-^ (r(t;a))' = 4, [T{t;d),P{t;d)] =0. (B.3) 



n=l 



From these operators, one can calculate the operator Q(t; C) (see Eq. fl2.14p ) which is 
given as 



The coefficients of the asymptotic expansion (12.181) are then calculated as 

(F-^C)^ 



(B.4) 



fiT-\ + 0(1/0, 



Y{t; = 4 + ^adj-^(F-2) [T-'{T-\ Hit)}] + 0(1/^), 



(B.5) 



where adj ^ is the inverse operator of adj(y4)[i?] = AB — BA. In the k = 3 case, by using 
the formula, adj"^(F-i)[X] = [T-\X]/3, one can show 



Y^{t)=^-(^H{t)-T-'H{t)T). 



(B.6) 



Here we have checked that ipo{t) = is true for first few cases k = 3,4,5 and this is 
consistent with our solutions. 

B.2 Practional-superstring {p,q) = (1,2) critical points (r = 3) 

In this case, we only study k = 2 case, but generally one can calculate as follows: The 
Lax operators in these cases are 

P{t; d) = rd + H{t), Qit; d) = {r-\t; d) P'{t; d)) ^ - fi{r-\t; d)) ^ 

= rd^ + H{t)d-S2{t)-fir-^ (B.7) 
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Therefore, the operator Q(t; () is given as 

Qit; = - r-'Hit) c - dH{t) - s^it) - ^t-\ 



or 



Q_3(t) = T-\ Q_2{t) = -T-'H{t), = ~dH{t) - S^it) - ^T-\ 

Here 52(t) satisfied 

[r, S2{t)] + TdH = 0, {r, ■ ■ ■ , r, S2{t)}k + {r, ■ ■ ■ , r, H{t), H{t)}k = o. 

In the k = 2 case, 5*2 (t) is given as 

1 



(B.8) 



(B.9) 



S2{t) = ^[cTif{t)-ta2df{t)y Hit) = la^fit). 
The coefficients of the asymptotic expansion are given as 

¥'(C) = ^i(y-K)+(^(l/C), 

y[C,) =12 + 7;^<^2— 7O-3— + -l(J2 7^ h U[l/Q ). 



(B.ll) 



(B.12) 



(B.13) 



C Supplements to Theorem [5] and Theorem [6 



In this appendix, we first show the explicit form of the linear expressions Eq. f l4.38p and 
Eq. f l4.48p . and then show some examples. Before we focus on these cases, we summarize 
the general properties of these systems by introducing the following four categories of the 
indices i of yn,i- 



(I) l<i< 



k + 3 
4 . 



(Ill) 'l±l + l<^< 



■.A, 

3k + 3 
4 . 



(II) B: = 
--: C, (IV) D : 



k + 3 



l<i< 



k + 1 



3k + 3 
4 . 



+ l<i<k, (C.l) 



which is closely related to the multi-cut boundary conditions Eq. f l4.33p and Eq. ( 14. 45 p . 
On the other hand, we show this division in the proffie i/^Y™'*- Here we show the categories 
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We define the symmetric product {^i, A2, ■ ■ ■ , Ak}k as 

k 



^ ^ {ill J 0-12 ; ■ ■ ' ''^jfclfe- 



(B.IO) 
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of I and III with italic font and the categories of II and IV with bold font: 
A; = 4A;o + l, (fcp g N) : 

^(sym) _ 



B 


(D 


D + 1) 


(A 


B + 1) 


{C 


D + 2) 


{A-1 


B + 2) 




(D 


B) 


(A 


D + 1) 


(c 


B + 1) 


{A-1 


D + 2) 


{C - 1 




D 


(A 


B) 


(C 


D + 1) 


{A- 1 


B + 1) 


{C - 1 


D + 2) 




. 


D) 


(C 


B) 


{A - 1 


D + 1) 


{C - 1 


B + 1) 


{A -2 







U 


2 


/fc+7 
V 2 


k) 


(5 


k+lN 
2 


; A 




{2 


k+3\ ■ 
2 1 




k-1) 


lk+1 
\ 2 


2 ^ 


(5 


k) 




k+lN 
2 ^ 




1) 


k+S 

9. 




\ 2 


k-1) 


{3 




Ik+b 
\ 9. 


k) 


{2 


khi) 


(k+3 

\ 9. 


1) 




k-3\ 
2 / 


(5 


k-1) 




kzl) 

2 / 




k) 


\ 2 


k+l\ 
2 / 


1 



and some concrete examples (A; = 9 and 13) are 



4 


(8 


9) 


{3 


5) 


(7 




{2 


6)^ 




(8 


4) 


{3 


9) 


(7 


5) 




1) 


6 




8 


{3 


4) 


(7 


9) 


{2 


5) 


{6 


1) 


■Jl 


[{3 


8) 


(7 


4) 




9) 


{6 


5) 


1 


■■Jo 



J2 
Jl 

Jo 

(C.2) 



<-'l3,2 



5 


(11 


12) 


(^ 


6) 


{10 


13) 


{3 


7) 


(5 


^) 


{2 


5)1 


(11 


5) 


(4 


12) 


(iO 


6) 


{3 


13) 


{9 


7) 


(^ 


1) 


5 


11 


(-^ 


5) 


{10 


12) 


(5 


6) 


(5 


13) 


{2 


7) 


(8 


^) 


L (4 
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5) 


{3 


12) 


(5 


6) 


(^ 


13) 


{8 


7) 


i 



J2 
Jl ■ 

Jo 

(C.3) 



k = 4:ko + 3, (feo e N) : 

^(sym) _ 
'-'k,2 — 



D 


(B 


B + 1) 


(C 


D + 1) 


(A 


B + 2) 


{C - 1 


D + 2) 




(B 


D) 


(C 


B + 1) 


(A 


D + 1) 


{C- 1 


B + 2) 


{A-1 




B 


(C 


D) 


(^ 


B + 1) 


{C- 1 


D + 1) 


[A - 1 


B + 2) 




AC 


B) 


(A 


D) 


(C- 1 


B + 1) 


{A - 1 


D + 1) 


{C -2 







(4 


k^) 


(¥ 


k) 


{3 




(k+5 
\ ^ 


i) 


{2 


fc+5\ - 




k-1) 


\ 2 


¥) 


(5 


k) 


V a 


k+l\ 
2 ) 


(^ 




fc+.'y 

■9 




(¥ 


k-1) 


(3 


— ) 

2 / 


(¥ 


k) 


(2 


k+l\ 
2 / 


(¥ 






k-3\ 
2 ) 


{3 


k-1) 




¥) 




k) 


V 2 


2 ^ 





and some concrete examples {k — 7 and 11) are 



7 


(3 


4) 


(^ 


1) 


{2 


^)1 


■Js 


(3 


7) 


(^ 


4) 


{2 


1) 


5 


■J2 


3 


(G 


7) 


(2 


4) 


(o 


1) 


■■Jl 


L(^ 
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7) 


(5 


4) 


1 


: Jo 



Jz 

J2 
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Jo 
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^(sym) 
>-^ll,2 



" 10 


(4 


5) 


(9 


11) 


(5 


6) 


{8 


1) 


{2 


7)1 


: ^3 


(4 


10) 


{9 


5) 


{3 


11) 


(5 


6) 


{2 


1) 


7 


: J 2 


4 


{9 


10) 


(5 


5) 


{8 


11) 
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: Ji 
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4) 
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(5 


5) 


{2 
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(7 
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i 


: Jo 



(C.5) 



If one follows trajectories of numbers, One may notice that the trajectories of the numbers 
in the region I and III (written by italic font) almost form slash shape, and the 

trajectories of the numbers in the region II and IV (written by bold font) almost form 
backslash shape, 



In the left-hand and right-hand ends of the profile J'lfl'^\ there 



are a few exceptions which form curved shape as "<" and so on. From this property, we 
can see the following facts: 

• Stokes multipliers si^ij O are almost given by i G II, IV (bold font) and 
j G I, III {italic font). We emphasize this fact by writing s^^ij ^ (j|i)/- Then there 
are only a few exceptions, s*,ij (i, j G II, IV) and (i, j G I, III), which appear in 
the left-hand and right-hand ends of the profile J'^l"^'' ■ Interestingly, there is no 
Stokes multipliers of the type s/ j j with i G I, III and j G II, IV in this r = 2 case. 

• From this fact, one can show that the relation between the symmetric Stokes mul- 
tipliers SqJ™^ and the fine Stokes multipliers si^ij (Eq. f l3.16p ) are almost trivial: 



,(sym) 



S3; j for almost all {j\i)i G Jkl'^\ and that only the following few multi- 



pliers are the exceptions: 



k + S 
2 



2' 



„(sym) 



Si k+i 

' 2 



g -|- Sq k+1 k±3_ Sg h±S_ g, (C.6) 



and 
A; = 4A;o + 1 : 



Jsym) 
'0,B.A 



= 4A;o + 3 : 

„{syin) 



(,sym; (sym) . q\ 

0,D,C - "5l,D,C + -32,0,3 So,B,C, "Sq^B+I.C " "S2,B+1,C + -33,3+1,3 So,B,C- 



• These relations are important not only because they are used in deriving the results 
in this appendix, but also because they provide a concrete example which guarantees 
the claim shown in Eq. fl3.37p . In particular, we expect that one can extend this 
discussion to the general systems of (A;, r; 7) which are controlled by the method 
developed in Section |3J 



C.l The explicit form of Eq. ( lirSSj ) and Eq. ( jinsl ) 



Below is the explicit form of the linear expressions Eq. fl4.38p and Eq. (14. 48 p . Note the 
function tik) appears in these formulas is given as 

1 (/t = 4A;o + 3) • ^^-^^ 
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The linear expression of Eq. ( 14.381) 

Region I: 

yn,i{{ym,l}rnez) = Vn+i-l,! ^ 0, i G (I), (C.IO) 

Region II: 



yn,B+j{{ym,l}meI.) = Vn+A+j,! + ^ 4,B+i-a,A-i+a+e(fc) ^ ^"+^-i-l+2a+e(fc),l + 

a=0 

i-1 

+ X] ^0J3+j-a,A-j+a+l+e{k) ^ 2/n+yl-i+2a+e(fc),l ) (CH) 



a=0 

Region III: 

yn,i{{yn.,l}mez) =0, i G (III). (C.12) 

Region IV: 

j 



yn,D+j{{ym,l}mez) = ^O^D+j-o-A-j+a ^ ?/n+A-j-l+2a,l 

a=0 

j-l+e(fc) 

+ ^ ^Oj:>+j-a,A-j+a+l ^ 2/n+A-i+2a,l- (C.13) 
a=0 

The linear expression of Eq. (14.481) 

Region I: 

yn,i{{ym,^^}m€z) =0, i G (I), (C.14) 

Region II: 

i 

yn,B+j{{y„,^h±3}mez) = '^S^^Jj^j^^^c_j^^ X y^^c-j-!^+2a,!^ + 

a=0 

+ X] ■^0,B+i-a,C-j+a+l ^ ^n+C-i+2a+l-*i2, ' (C.15) 
a=0 

Region III: 

yn,i{{y„,^h^]mel) = y„+i-i±3,fe±3 7^ 0, i G (III), (C.16) 

Region IV: 



yn,D+j{{y^^k±l}rnez) = y„+D+j-±±2 ,±±3 + 



i-1 



I \ ^ (sym) — . 

+ / .^O.D+i-a.C-i+a+l-dk) ^ ?/n+C-j+a+l-e(fc)-i±2, *!±3 + 



a=0 



+ ■^0,D+j-a,C-j+a+2-e(k) ^ yn+C-j+a+2~e(k)-i^ 



2 ' 2 

a=0 



(C.17) 
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C.2 Some examples of Theorem [5] 

Below we show the concrete expressions of the equations in Theorem |5] for some special 
cases {k = 5, 7, 9 and 11). First of all, the vectors are expressed only by using {yn,i}n&z- 

( yn,i \ 

Vn+l,! 



( yn,i \ 



k = 5: 



■^0,4,2 yn,l 







k = 7: F^") 



, (sym) . 
Vo,5,2 Vn+l,!/ 



( 



yn+2,1 
(sym) 
'■^0,5,2 VnA 








k = 9 



yn,i 

yn+1,1 
yn+2,1 



\ ■'0,7,2 "T "^0,7,3 yn+2,1/ 

\ 



I (sym) 
l/n+3,1 + -50,4,3 l/n+2,1 
(sym) 
"■^0,6,2 yn,l 








■^0,8,3 yn+2,1 
\'5o,8,3 l/n+3,1 + Sq q^^ yn+2,1 + ^0,9,2 l/n+l,l/ 



A; = 11 : 



y(n) 



l/n,l 
yn+1,1 

yn+2,1 



\ 



yn+3,1 
I (sym) I (sym) 

l/n+4,1 + 50,5,3 2/n+2,l + "^0,5,4 yn+3,1 
(sym) 
~"^0,7,2 yn,l 









,(sym) 



,(sym) 



■'CICS 2/n+2,l + -50^10^4 yn+3,1 
\ „(sym) I „(sym) I (sym) (sym) / 

\*0,10,3 yn+3,1 "T ■'0,10,4 yn+4,1 "T ■^0,11, 2 yn+1,1 "T ■^0,11,3 yn+2,l/ 



(C.l^ 



Secondly, the multi-cut BC recursions are expressed as 
k = 5 



'^5[yn,l] 
^5[yn,l] 

k = 7 

J^7[yn,l] 

Qiiyn,^ 

k = 9 

'^9[yn,l] 
^9[yn,l] 



■ yn+2,1 + -51,3,2 yn+1,1 + -53,4,2 2/n,l — 0, 
-S0,5,2 l/n+2,1 + -52,1,2 l/n+1,1 " yn,l = 0, 

■ yn+3,1 + -53,4,3 yn+2,1 + -5l,4,2 2/n+l,l + -53,5,2 2/n,l = 0, 
-52,7,3 l/n+3,1 + -5o,7,2 l/n+2,1 + -52,1,2 yn+1,1 - yn,l = 0, 

yn+4,1 + -51,4,3 yn+3,1 + -53,5,3 yn+2,1 + -5l,5,2 yn+1,1 + -53,6,2 yn,l = 0, 
-50,8,3 yn+4,1 + -52,9,3 yn+3,1 + -5o,9,2 yn+2,1 + -52,1,2 yn+1,1 "~ yn,l = 0, 
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k = ll 

^ n[yri,l\ = l/n+5,1 + -53,5,4 ?/„+4,l + Si,5,3 7/„+3,i + 83,5,3 l/n+2,1 + ■5l,6,2 ?/n+l,l + -53,7,2 ?/n,l = 0, 

= -52,10,4 Z/n.+5,l + -5o,10,3 Z/n+4,1 + -52,11.3 l/n+3,1 + -5o.ll,2 Z/n.+2,l + -52,1,2 Z/n.+ l,l "~ l/n,l = 0. 

(C.19) 

D Derivation of the continuum solutions 

In this subsection, we derive the continuum solutions of Theorem [HI According to Lemma 
[H the monodromy free condition can be solved if the matrix Sq^^"^'' F"^ is diagonalizable. 
For the continuum solutions, we solve this problem by requiring that all the eigenval- 
ues of the matrix Sq'^^^ are distinct. This means that we require the characteristic 
polynomial Tilx) of the matrix S^^^^ F~^ satisfy 

^(sym) p-l\ _ k 



H{x) = det [xlk - T-'j =x' -1. (D.l) 

The coefficients of the characteristic polynomial are related to the Stokes multipliers and 
here are several examples: 

The 5- cut case: 

= -1 + + X (-Si,3,2 + -50,5,2^2,3,5 " ^3,1,5) + 
+ (-50,3,4-50,5,2 ~ -5l,l,4 + -52,1,2-52,3,5 ~ -51,3,2-53,1,5 ~ -53,4,2) + 
+ X^ (--50,5,2 - -51,1,4-51,3,2 + -50,3,4-52,1,2 - -52,3,5 - -53,1,5-53,4,2) + 
+ ( — -50,3,4 — -52,1,2 ~ -51,1,4-53,4,2) ; 

The 7-cut case: 

n{x) = -l + x"^ + X (-Si,7,6 + -50,3,6-52,7,3 " -53,4,3) + 

+ (-50,3,6-50,7,2 — -5l,4,2 + -52,4,6-52,7,3 ~ -53,1,6 ~ -51,7,6-53,4,3) + 

+ X^ ( — -51,1,5 ~ -51,4,2-51,7,6 + -50,3,6-52,1,2 + -50,7,2-52,4,6 + -50,4,5-52,7,3 ~ -53,1,6-53,4,3 — -53,5,2) + 
+ X^ ( — -50,3,6 + -50,4,5-50,7,2 + -52,1,2-52,4,6 " -52,7,3 — -51,4,2-53,1,6 — -51,1,5-53,4,3 " -51,7,6-53,5,2) + 
+ X^ ( — -5o,7,2 ~ -51,1,5-51,4,2 + -50,4,5-52,1,2 ~ -52,4,6 ~ -53,1,6-53,5,2) + 

+ X^ (--50,4,5 - -52,1,2 - -51,1,5-53,5,2) • (E).2) 

These equations become simpler if one uses the notation given in Eqs. (14.571) and fl4.58p . 
One can read the general formula for k = 2m + 1: 

m n n 

n=l i=l 1=0 

m n n 

+ E [E ^i^n-i - E ^m+l-iOm^n+i , (E).3) 

n=l 1=0 1=1 
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where we have introduced 6q = 1. Therefore, by comparing both sides of Eq. (ID.ip . we 
obtain constraints on the Stokes multiphers: 

= ^m^m -01- 6i, 

n 0*0* I o* o* o o 

n 0*0* \ o* o* \ o* o* o o o o 

n 0*0* \ o* o* \ o* o* \ o* o* o o o o o o o o 



Since a half of the Stokes multiphers {6'„}^^i are given as 



(D.4) 



(D.5) 



we fix all the other Stokes multipliers from these constraints. Note that, since 

all the eigenvalues are distinct, the indices (ni,n2, ■ ■ • ,nm) for {6'„} are also m (= [|J) 
distinct integers. Here we can freely choose the ordering: 

(rii, 712, • • • , : \ <nx < < ■ ■ ■ < <k. 

With noting the following relation: 

k 

.2. 



0*0 o* 



m 



(D.6) 



(D.7) 



and recursively rewriting the constraint for the continuum solutions, Eqs. ( ID.ip and (ID.Sp . 
we obtain the following simple form: 

On = S4{9,},^^) + e*^_^^,e*^, (n = 1, 2, ■ ■ ■ , m) , (D.8) 

with the polynomials Sn{x) (defined by Eq. (14.811) ). This results in Theorem [81 



E Calculation in the 3-cut (1, 1) critical point (r = 2) 

The specialty of the 3-cut (1, 1) critical point is that the symmetric Stokes sectors Z)4„ 
(see Eq. (13.151) ) do not cover the whole plane C. Therefore, we consider doubling of the 
sectors 



2ny 



Q(sym) 
^2n 



= 'S'2„S'2„+l, 



[n 



0,1, 



5), 



and express the boundary condition (14.331) as follows: 

y(4n) ^ 



y4n,l 



(4n) 



X 



n+2 



V 



y(4ri+2) 



^l/4n,3/ 
yin+2,1 ' 

Z/4n+2,2 I = r"X(^"+2) 
,Z/4n+2,3> 



/ ^(4n+2) 
V ^(4"+2) 



0, 



(E.l) 



(E.2) 
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with 



^(2n) ^ ^(«^y-)x(2«+2), (n = 0, 1, ■ ■ ■ , 5). 

This is then written as 

y(4n) _ ^^sym)y(4„+2) y(4n+2) _ ^^(sym)p_l-j y (4„+4) 

y4n,i — y4n+2,i + Z^j Sq ^ j X ?/4„+2,j , 

I v-^ (sym) ] 

yin+2,i — yAn+A,i-l + l^j S2,j,j ^ yAn+A,j-l\ ■ 

These recursion relations are expressed as 

l/4n,3 = yAn+2,3 = 0, yAn,l = Z/4n+2,l 7^ 0, 

yAn,2 = yAn+2,2 Z/4n+2,2 — yAn+A,l 

and the following two recursion equation for y4n,i 

2/4n,l = ■52,1,2 X ?/4n+4,l, l/4n+4,l = —-53,3,2 X l/An,!- 



(E.3) 



(E.4) 



As one may notice, this equation itself is the same as Eq. (14.371) . The solutions 
by /) to this boundary condition is easily solved as 



JO 

'3,3,2 



JO 

'2,1,2 



(/ = 0,1,2), 



and the general solution is given as 

■'0,2,3 



(E.5) 

(E.6) 
(labeled 

(E.7) 
(E.8) 



with Eq. ( ]E.7p . This provides the first case of the continuum solution ( ID.SP . 



F Calculation in the 4-cut (1, 1) critical point (r = 2) 

Here we calculate the 4-cut (1, 1) critical point as an example in which the coprime 
condition of Eq. (13. 8p is violated: 



(fc,r) = (4,2). 

In this case, the leading exponents are degenerate: 
and we consider the subleading Stokes lines: 



Re 



f (1) (3) \^r-l 



(3) 



0, 



Re 



f (2) (4j \/-r-l 



(4) 



0. 



(F.l) 



(F.2) 



(F.3) 



The dominance profile in the ( plane is shown in Fig. [121 

Here we use the fine Stokes sectors D„ (calculated in the leading Stokes lines) which 
are defined as 



n = 0, 1,2,3. 



(F.4) 
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Figure 12: The dominance profile in the 4-cut (1, 1) case in terms of (. The bold lines express the 
leading Stokes lines with degeneracy ip^^'^ ^ ip^'^'^ and The dashed lines express the sub 

leading Stokes lines for (1,3) and (2,4). 



All fine Stokes matrices can be expressed in terms of Sq as 



Sn = r-"5or", 



So 



/I \ 

a 1 /3 

e 1 

V7 5 1/ 



(F.5) 



(F.6) 



Then the multi-cut boundary condition is given as 

(yn,i ^ o\ 

y(n) _ ■p?ij^(n) Z/re,2 = 

~ yn,z = 

\yn,4 7^0/ 

The recursive equations are expressed as 



yn,l = yn+lA, = e X y„+i,4, yn+1,1 + a X = 0, 7 X Vn+i^i - ?/„,i = 

(F.7) 



and there are four solutions which are labeled by / (a — )■ a 



(0^ 



e« = 0, y£\ = c"' 



(/ = 0,1,2,3). (Fi 



By directly solving the monodromy free condition, the other Stokes multipliers (Z?*^'-' and 
7*^')) are also fixed and the solution is given as 



S() 



1 1 




\ 


-w' 1 


-co-' 







1 








V 



(/ = 0,1,2,3). 



(F.9) 
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